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ed  at  Project  MAC,  MIT.  A  typical  problem  in  applied  mathematics,  namely  asymptotic 
analysis  of  a  definite  integral,  is  solved  using  MACSYMA  to  demonstrate  the  usefulness 
of  such  a  system  and  the  facilities  provided  by  WANDERER. 
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This  result  gives  the  first  four  terms  of  the  desired 
asymptotic  expansion  in  C. 
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EVALUATION  OF  DEFINITE  INTEGRALS 
BY  SYMBOLIC  MANIPULATION* 


Abstract 


A  heuristic  computer  program  for  the  evaluation  of  real 
definite  integrals  of  elementary  functions  is  described.  This 
program,  called  WANDERER  (WANg's  DEfinite  integRal  Eva'iuatoR), 
evaluates  many  proujr  and  improper  integrals.  The  improper 
integrals  may  have  a  finite  or  infinite  range  of  integration. 
Evaluation  by  contour  integration  and  residue  theory  is  among 
the  methods  used.  A  program  called  DELIMITER  (DEfinitlve  LIMIT 
EvaluatoR)  is  used  for  the  limit  computations  needed  in  evaluat¬ 
ing  some  definite  integrals.  DELIMITER  is  a  heuristic  program 
written  for  computing  limits  of  real  or  complex  analytic  func¬ 
tions.  For  real  functions  of  a  real  variable,  one-sided  as  well 
as  two-sided  limits  can  be  computed.  WANDERER  and  DELIMITER  have 
been  implemented  in  the  MACSYMA  system,  a  symbolic  and  algebraic 
manipulation  system  being  developed  at  Project  MAC,  MIT.  A  typi¬ 
cal  problem  in  applied  mathematics,  namely  asymptotic  analysis  of 
a  definite  integral,  is  solved  using  MACSYMA  to  demonstrate  the 
usefulness  of  such  a  system  and  the  facilities  provided  by 
WANDERER. 
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CHAPTER  I 
INTRODUCTION 

In  recent  years  rapid  advancement  has  taken  place  In 
the  art  of  using  modern  electronic  computers  to  facilitate 
symbolic  mathematical  computations.  Computer  based  software 
systems  for  this  purpose  are  generally  known  as  symbolic 
manipulation  systems.  Such  a  system  Is  designed  to  aid  Its 
users  In  handling  mathematical  expressions  and  functions. 

As  the  field  of  symbolic  and  algebraic  manipulation  grows 
stronger  and  more  sophisticated  so  do  the  various  algebraic 
manipulation  systems.  MACSYMA  [16]  Is  a  recent  redesign  of 
earlier  systems  Incorporating  many  new  Ideas  and  results 
developed  In  the  field.  The  development  of  MACSYMA  began  at 
Project  MAC/  MIT/  In  early  1969,  The  effort  has  been  guided 
by  Professors  Martin  and  Moses  of  MIT  and  Involves  a  group 
of  researchers  and  students  Including  the  author.  Drawing 
on  the  past  work  of  Martin  [15]/  Moses  [20]  and  Engelman 
[9]/  MACSYMA  soon  evolved  Into  a  system  which  extended  the 
capabilities  of  automated  algebraic  manipulation  systems  In 
many  new  areas.  Among  these  are  the  ability  to  compute 
limits  of  functions  and  the  ability  to  evaluate  definite 
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Integrals. 

The  limit  concept  Is  basic  to  mathematical  analysts. 
Being  able  to  compute  limits  automatically  greatly  Increases 
the  potential  of  a  symbol  manipulation  system  In  doing 
analytical  mathematics.  In  fact  the  evaluation  of  definite 
Integrals  Is  heavily  dependent  on  the  limit  process,  as  Is 
the  expansion  of  functions  In  power  series  and  many  other 
mathematical  problems.  Although  the  computation  of  limits 
has  been  studied  previously  to  some  extent  [11,  13],  we 
describe  a  limit  program  called  DELIMITER  (DEflnltlve  LIMIT 
EvaluatoR)  which  Is  more  powerful  than  previous  programs. 

It  Is  discussed  In  detail  In  chapter  2, 

The  problem  of  computing  Indefinite  Integrals 
symbolically  by  computers  has  been  Investigated  rather 
thoroughly.  First  among  the  computer  programs  developed  for 
this  purpose  was  SAINT  (Symbolic  Automatic  INTegrator)  by 
Slagle  In  1961.  A  more  powerful  program  named  SIN 
(Symbolic  INtegrator)  [23]  was  developed  In  1967  by  Moses. 
Theoretical  work  In  this  area  Include  Richardson's 
undectdabl 1 1 ty  result  for  a  certain  classes  of  Integrals 
[22]  and  Rlsch's  decision  procedure  for  determining  the 
existence  of  the  Indefinite  Integral  of  a  member  In  the 
class  of  elementary  functions  [23,  24,  25,  26],  There  Is  a 
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comprehensive  review  by  Moses  on  the  progress  during  the 
past  ten  years  In  this  particular  area  of  symbol  manipula¬ 
tion  [21]. 

However/  hardly  any  work  has  been  done  In  the  direction 
of  definite  Integration.  Reported  here  Is  a  first  definite 
Integration  program  called  WANDERER  (WANg's  DEflnlte 
IntegRal  EvaluatoR).  Experiments  with  DELIMITER  and 
WANDERER  are  the  principal  subjects  of  this  thesis. 

WANDERER  has  been  Implemented  In  MACSYMA  and  makes  use  of 
many  facilities  provided  In  It.  Some  of  these  facilities 
are:  Input  and  two-dimensional  output/  slmpl If tcatlon, 

solution  of  polynomial  and  systems  of  linear  equations/ 
canonical  rational  function  simplification  and  the  SIN  [20] 
Integration  program. 

Evaluation  of  definite  Integrals  can  sometimes  be  as 
easy  as  computing  the  Indefinite  Integral  of  the  given 
Integrand  then  substituting  In  the  limits  of  Integration. 
WANDERER  computes  Integrals  of  this  kind  by  using  those 
parts  of  SIN  that  have  being  Implemented  In  MACSYMA. 

However/  many  Interesting  definite  Integrals  are  not 
obtainable  In  this  manner.  In  some  cases  they  may  be 
Improper  Integrals  or  their  Indefinite  Integrals  do  not 
exist.  In  other  cases  It  Is  easier  to  evaluate  the  definite 


Integral  directly  than  to  obtain  the  Indefinite  Integral 
f I rst. 
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A  few  limit  and  Integration  problems  are  listed  below 
to  give  an  Indication  as  to  the  scope  of  WANDERER  and 
DELIMITER.  The  results  obtained  by  the  programs  can  be 
found  In  chapters  2,  4  and  5. 


roo  2  cq  6 

fx+AX+B  f  XdX 

/  . <1X  /  . 

/  4  2  L  2  15/2 
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integrals  with  an  Infinite  range  are  discussed  tn 
chapter  4,  those  with  a  finite  range  In  chapter  5. 
Important  routines  and  algorithms  are  collected  In  chapter 
6. 


An  application  of  the  MACSYMA  system  to  an  asymptotic 
analysis  problem  Is  Illustrated  In  chapter  7.  The  purpose 
is  to  show  how  such  a  symbol  manipulation  system  can  be  used 
to  solve  non-tr!vlai  problems  that  may  occur  frequently  In 
applied  mathematics.  This  application  demonstrates  the 
usefulness  of  the  many  faculties  provided  by  MACSYMA  and 
WANDERER. 

A  timing  experiment  has  been  conducted  to  check  the 
performance  of  WANDERER  and  DELIMITER.  The  results  are 
Included  In  appendix  C.  Because  of  the  limited  character 
set  of  most  computer  consol es,  some  special  symbols  are 
needed  to  denote  the  frequently  used  mathematical  functions, 
constants  and  operators.  For  example,  *  and  **  are  used  to 
denote  multiplication  and  exponentiation  respectively. 
Appendix  D  contains  a  list  of  notations  we  shall  use. 

From  Richardson's  undecldabll Ity  results  [22],  we  have 
shown  that  the  convergence  of  a  class  of  Integrals  of 
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elementary  functions  Is  recursively  undecldable.  The  proof 
ts  In  appendix  E. 

The  remainder  of  this  chapter  serves  as  an  Introduction 
to  algebraic  manipulation  systems.  If  the  reader  is 
familiar  with  such  systems  he  may  proceed  directly  to  the 
next  chapter. 

An  algebraic  manipulation  system  Is  a  collection  of 
computer  programs  designed  to  facilitate  the  solutton  of 
mathematical  problems.  Such  a  system  has  the  ability  to 
handle  both  symbols  and  numbers.  This  capability  to  mani¬ 
pulate  symbols  mathematically  Is  what  differentiates  these 
systems  from  the  various  computer  subroutines  which 
specialize  In  numerical  analysts. 

In  using  symbol  manipulation  systems  such  as  these,  one 
usually  Interacts  with  It  In  a  time-sharing  environment  via 
a  typewr I ter-1 Ike  console.  Normally,  data  and  commands  are 
typed  In  by  the  user.  Results  from  the  computer  are  sent 
back  to  the  console  for  display.  This  arrangement  attempts 
to  provide  a  user  the  ease  and  flexibility  of  the  pencil  and 
paper  he  Is  so  accustomed  to,  whtle  permitting  the  computer 
to  assist  him  In  his  algebra  and  formal  deduction  from  one 
step  to  the  next.  If  the  manipulations  Involved  are  non- 
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trivial,  so  much  the  better.  Some  of  the  manipulations  such 
a  system  can  provide  are:  GCD  calculations,  factoring,  ra¬ 
tional  function  arithmetic,  matrix  manipulation,  solution  of 
algebraic  equations,  solution  of  systems  of  linear  equa¬ 
tions,  Indefinite  Integration  of  elementary  functions* 

Let  us  take  a  closer  look  at  MACSYMA  as  a 
representative  model  of  other  systems.  MACSYMA  receives 
Inputs  In  the  form  of  linear  character  strings  typed  by  the 
user.  FORTRAN-1  Ike  notation  Is  used  for  the  Input.  For 
example. 

The  user  types:  4*(X**2+X+1)/Y@ 

The  @  sign  signifies  the  end  of  a  command  string.  As  a 
result  of  this  command,  a  two  dimensional  display  Is 
returned. 

2* 

4  (X  +  X  ♦  1) 

MACSYMA  types:  . - 

Y 

Inside  MACSYMA,  expressions  are  represented  by  list 
structures  In  a  prefix  notation  (common  to  many  systems). 

For  example  the  expression  above  would  be  represented 
Internally  by  something  like 

(TIMES  4  (PLUS  1  X  (EXPT  X  2))  (EXPT  Y  -1>) 


Inside  MACSYMA 
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Other  forms  of  Internal  representations  are  possible. 
The  choice  of  Internal  representation  Is  an  Important  aspect 
of  the  design  of  a  system.  Of  course  more  than  one  form  of 
Internal  representation  can  be  employed  In  a  single 
algebraic  manipulation  system.  In  fact  MACSYMA  has  a 
special  Internal  representation  for  rational  functions  which 
Is  used  to  gain  efficiency  during  certain  polynomial  and  ra¬ 
tional  function  manipulations. 

One  simple  application  of  MACSYMA  Is  factoring  the 
polynomial 

3  2 

P(X)  *  X  +  U  -11  X  -  30 

A  user  who  wants  to  factor  PCX)  using  MACSYMA  types: 

FACTORCX**3+4*X**2-ll*X-33)@ 

This  Input  command  causes  the  factorization  of  PCX)  over  the 
Integers.  The  output  Is 

(X  +  2)  (X  -  3)  (X  ♦  5) 

While  this  problem  may  seem  easy#  factorization  of 
polynomials  of  higher  degree  can  be  very  difficult  to  do  by 
hand.  Indeed/  algebraic  manipulation  systems  can  be  most 
helpful  when  one  wants  to  manipulate  complicated  functions 
and  expressions.  MACSYMA  can  carry  out  accurately  with 
great  speed:  summation  over  IndtceS/  expansion  of  products 
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and  powers/  calculation  of  large  determinants/  Inversion  of 
large  matrices. 

Applications  In  pure  mathematics  Include  computations 
In:  number  theory/  group  theory/  Lie  algebras  and  set 
theory. 

Research  In  the  field  of  symbol  and  algebratc  manipula¬ 
tion  has  led  to  many  new  results.  The  fast  GCD  algorithm 
[4]  and  the  finite  field  arithmetic  polynomial  factoring 
algorithm  [2]  are  two  examples.  Continued  work  In  this  area 
will/  hopefully/  result  In  computer  systems  which  are 
Increasingly  valuable  to  engineers  and  mathematicians. 
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CHAPTER  tl 
DELIMITER 


0.  Introduction 

The  limit  concept  Is  fundamental  to  mathematical  analy¬ 
sis.  Baste  concepts  such  as  the  rules  of  differentiation 
are  derived  from  limiting  processes.  More  complex  problems 
such  as  Improper  Integration,  convergence  of  series,  series 
expansion  of  functions  and  contour  Integrals,  to  name  a  few, 
also  require  the  computation  of  limits  In  their  solution 
process. 

Therefore,  one  can  expect  that  automatic  computation  of 
limits  would  greatly  Increase  the  capability  of  a  symbolic 
mathematics  system  In  doing  analytical  mathematics.  The 
programs  described  below  provide  such  a  capability  In  the 
MACSYMA  system  [16]. 

Automating  the  computation  of  limits  has  been  studied 
previously  to  some  extent.  Fenlchel  [11]  discussed  certain 
decidability  problems  of  limits  and  provided.  In  the  FAMOUS 
system,  some  basic  routines  for  computing  two-sided  limits. 
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In  his  thesis,  Iturrlaga  [13]  worked  on  both  one  and  two- 
sided  limits.  In  addition  to  (.'Hospital  rule,  he  applied 
asymptotic  analysts  to  quotients  of  polynomials  essentially 
by  replacing  those  polynomials  by  their  leading  terms.  This 
can  be  done  when  the  variable  In  the  limit  approaches  Infi¬ 
nity  so  that  the  leading  term  Is  asymptotically  equivalent 
to  the  original  polynomial.  Moreover,  he  discussed  limits 
of  sequences. 

DELIMITER  Is  more  powerful  than  either  Fenlchel's  or 
Iturrlaga's  limit  program,  L'HospItal's  rule  Is  a  basic 
method  used  for  Indeterminate  forms.  In  addition,  this  pro¬ 
gram  employs  a  fast  routine  for  limits  of  rational  func¬ 
tions.  It  has  an  efficient  algorithm  for  a  class  of  expres¬ 
sions  called  RP-express Ions  which  Involve  radicals  of  poly¬ 
nomials.  It  also  applies  the  method  of  reducing  complicated 
expressions  by  replacing  subexpressions  with  asymptotically 
equivalent  expressions.  The  method  of  comparing  orders  of 
Infinity  of  expressions  and  several  other  heuristic  methods 
are  used.  In  some  cases,  power  series  expansions  are  employ¬ 
ed  to  obtain  the  limits. 

This  chapter  Is  based  on  a  paper  by  the  author  [29] 
presented  at  The  Second  Symposium  on  Symbolic  and  Algebraic 
Manipulation,  Los  Angeles,  March,  1971. 


,  i«*  t-  -''r' 1  •'U-  $".  ;  '4*  )'*•■  T  'f '  W'  \  •  ,*'■'  '• 
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1.  Definitions  and  Symbols 

DELIMITER  Is  designed  for  finding  the  limit  of  a 
single-valued  function  f(X)  of  a  real  or  complex  variable  X 
as  X  approaches  some  limit  point.  One-sided  limits  can  be 
computed  If  X  and  f (X )  are  real-valued.  The  classes  of 
functions  allowed  by  this  program  Include  rational,  radical, 
logarithmic  and  exponential  functions,  and  also  the  trigono¬ 
metric  functions  SIN,  COS  and  TAN,  and  the  hyperbolic 
functions  SINH,  COSH  and  TANH.  The  general  form  of  a  user 
command  Is 

L I M I T ( ex  p  re  ss I on , variable, value,  dlrectl on ) @ 

with  the  fourth  argument  optional.  PLUS  as  a  fourth  argu¬ 
ment  Indicates  the  one-sided  limit  from  above,  MINUS  from 
below.  The  absence  of  the  fourth  argument  Indicates  no 
restriction  on  the  direction  of  approach.  One-sided  limits 
are  not  allowed  for  complex-valued  functions  or  limit 
points. 

If  f(X)  Is  not  continuous  at  X=a,  the  two-sided  limit, 
LIMlTCf (X),X,a),  does  not  exist.  In  such  a  case  the  symbol 
UNO  Is  returned  by  the  program  as  an  answer. 
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The  following  are  some  symbols  used  by  this  program. 


ZER0+ 

0  + 

ZERO- 

0- 

UND 

undefined 

IND 

Indefinite 

but  bounded 

INFINITY 

the  point 

at  Infinity  In 

the  complex  plane 

INF  and  MINF  wMl  continue  to  Indicate  respectively 
positive  and  negative  Infinity.  Here  are  some  examples  of 
how  these  symbols  are  used  Internally. 


LIMIT(COS(X),X,*Pl)  «  -1+  ZERO* 

LIMIT(SIN(X),X, INF)  »  IND 

L I M I  T( 1/S  I N(X),X, INF)  =  UND 

L I M I T  C 1/ ( X  —  %  I ) ,  X  ,  $  I )  -  INFINITY 
The  direction  from  which  a  limit  point  Is  approached  Is 
Important.  By  use  of  ZER0+  and  ZERO-,  simplification  rules 
such  as 

l/(-l  ♦  ZER0+  ♦  1)  — ►  INF 

and 

l/(-l  ♦  ZERO-  ♦  1)  — >  Ml NF 
are  possible  In  DELIMITER. 
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2.  Baste  Rules  and  Outline  of  Algorithm 

Some  of  the  rules  used  In  this  program  are  the  so-call¬ 
ed  "trivial"  ones  for  limits  of  continuous  functions# 
namely#  the  limit  of  a  sum  Is  the  sum  of  the  limits#  etc. 
There  Is  a  limit  routine  for  each  of  the  functions  SIN#  COS# 
TAN#  SINH#  COSH#  TANH  and  LOG  (base  %l) .  For  other  func¬ 
tions  the  rule 

LI  Ml  TC  f (g(X))#X#L>  «  LIMIT(f(Y)#Y#LIMIT(g(X)#X#L)) 

Is  used. 

For  one-sided  limits  where  the  variable  approaches  a 
point  other  than  0#  a  change  of  variable  Is  made  to  bring 
the  point  to  the  origin.  For  example#  the  limit 

LIMIT(f(X)#X#a# PLUS ) 

Is  converted  to 

LIMlT(f (a+Y)#Y#ZER0+) 

There  Is#  of  course#  a  complete  set  of  rules  governing 
simplification  of  the  new  symbols.  They  are  all  of  the  fol¬ 
lowing  nature 
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IND/INF  »  0 
-  ZERO+  =  ZERO 
1/ INFINITY  »  0 
Ml NF**2  =  INF 
LOG( C+ZERO+ )  » 


-1/INF  »  ZERO- 
-  INF  »  MINF 
-1/INFINITY  -  0 
(ZERO+)+ZERO-  *  0 
LOG(C)  ♦  ZERO* 


lie  will  now  briefly  outline  the  algorithm/  and  then  go 
Into  the  details  of  some  of  the  component  routines.  Let  us 
cons Ider 


LIMIT  (E(X)/X/L) 

where  L  Is  any  number  or  symbol  Including  INFINITY/  INF/ 
ZER0+  and  ZERO-.  If  L  “  MINF/  It  Is  set  to  INF  by  a  change 
of  variable  Y=-X.  Upon  receiving  the  arguments/  namely 
E( X) /  X  and  l,  DELIMITER  checks  whether  L= INFINITY  or  E(X) 
Involves  %!.  If  either  or  both  Is  true  then  a  global 
Indicator  CPLX  Is  set  to  the  value  TRUE  which  Indicates  that 
the  given  limit  Is  to  be  evaluated  over  the  field  of  complex 
numbers.  Otherwise  CPLX  Is  set  to  the  value  FALSE  which 
signifies  a  limit  problem  over  the  reals.  If  CPLX=TRUE/  all 
use  of  INF  and  MINF  are  replaced  by  INFINITY  and  the  notion 
of  approaching  a  limit  point  from  one  side  Is  no  longer 


valid 


‘  t  *T-  r*  ■ 
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The  following  diagram  Is  a  simplified  overview  of 
DELIMITER. 


DELIMITER 

SIMPLE  LIMIT - ►RETURN  ANSWER 


RATIONAL  FUNCTION - ^RETURN  ANSWER 


RP-FORMS- - ^RETURN  ANSWER 


BASIC  LIMIT  L' HOSPITAL'S  SERIES  ORDER 

RULES  RULE  EXPANSION  OF  INFINITY 


Fig.  I 

As  Indicated  In  fig.  1,  DELIMITER  has  four  stages. 
Results  of  simple  limit  problems  are  Immediately  returned  by 
the  first  stage.  Simple  limits  Include  two  cases:  1)  E  Is 
Independent  of  X.  2)  E(X)  =  X.  In  the  second  stage#  a  fast 
algorithm  Is  used  to  obtain  limits  of  rational  functions. 

The  third  stage  employs  an  algorithm  for  limits  of  RP-forms. 
If  the  first  three  stager;  do  not  produce  an  answer#  the 
given  limit  problem  enters  the  fourth  stage  which  contains  a 
variety  of  methods  Including  the  four  principal  ones  shown 
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In  the  figure.  In  this  stage  the  program  will  first  try  to 
extract  the  numerator  N(X)  and  the  denominator  D(X)  of  E(X) 
so  that  E ( X )  «  N(X)/0(X).  If  both  U(X)  and  0(X)  are  RP- 
forms  the  limit  Is  computed  by  an  algorithm  specially  for 
these  forms  (see  section  5  for  definition  and  algorithm). 
Otherwise  the  rule 

1 1m  E ( X )  »  1  Im  N(X)  /l Im  D(X) 

Is  applied.  If  both  Urn  N ( X )  and  11m  D(X)  are  0  or  Infinite 
we  have  an  Indeterminate  form.  In  this  case  L'HospItal's 
rule  will  be  applied  to  E(X)  with  one  exception:  when  L  Is 
INF  and  11m  N( X )  and  llm  D(X)  are  both  Infinite  and  both 
N(X)  and  D(X)  contain  exponential  functions  of  X  which  tend 
to  INF  as  X  approaches  L.  In  this  case  the  method  of  com¬ 
paring  orders  of  Infinity,  to  be  described  shortly.  Is  ap¬ 
plied.  If  DELIMITER  can  not  find  an  answer  Indication  of 
failure  will  be  returned. 

Sometimes  the  program  needs  to  know  the  value  range  of 
a  symbolic  parameter  In  order  to  compute  the  limit.  In  such 
cases,  the  program  will  query  the  user  at  his  console. 

It  Is  Important  to  note  that  these  programs  do  not 
store  any  table  of  limits  Thus  every  limit  obtained  Is  a 
result  of  computation. 
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Other  methods  available  Include  :  change  of  variable, 

discontinuity  tests  and  analysts  of  the  behavior  of  a  func¬ 
tion  near  a  finite  point.  A  flowchart  In  appendix  A  pre¬ 
sents  the  flow  of  control  In  a  more  complete  manner. 


^ - ^r**nv»V*w  *•>*• «?-*•■- r*  --'■"M  -*wv  >. ' •  :,  •***-»**/ eM  \-»ww1 
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3.  The  Limit  of  a  Rational  Function 

Let  P(X)  and  Q(X)  be  polynomials  In  X  and  possibly 
other  variables  and 

E( X)  =P(X)/Q(X) , 

n  »degree(P)-degree(Q), 

c  "leading  coeff lc!ent(P)/leadlng  coeff lclent(Q), 

Consider  LI  Ml T( E(X) /X/ L)  where  E(X)  or  X  may  be  complex 
valued.  The  following  algorithm  Is  used  to  compute  the 
1 Iml t. 

(1)  If  L  Is  finite  then. 

If  P(X)  *  Q(X)  3  0  then  set  P  and  Q  according  to 
P(X)»  P(X>  /(X-L),  Q(X)*  QCX)  /(X-L), 
by  long  division,  then  go  to  step  (1); 

If  Q(L)  Is  not  0  then  the  answer  Is  E(L); 
otherwise,  the  answer  Is  INFINITY  If  CPLX=TRUE 
and.  If  CPLX  Is  FALSE,  the  result  Is  INF  or  MINF 
depending  on  the  direction  from  which  Q(X)  ap¬ 
proaches  0  as  X  approaches  L,  and  of  course  the 
sign  of  P(L). 


(2)  If  L  Is  Infinite, 
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If  n=0,  the  answer  Is  c; 

If  n  <  0,  then  If  CLPX  =  TRUE,  answer  Is  0,  but  If 
CPLX  Is  FALSE  then  the  answer  Is  ZERO+  provided 
the  sign  of  c*L**n  Is  ♦  and  ZERO-  In  case  the  sign 
of  c*L**n  Is 

otherwise,  (n  >  0)  the  answer  Is  INFINITY  If 
CPLX=TRUE,  and  the  answer  Is  INF  or  MINF, 
depending  on  the  sign  of  c*L**n,  If  CPLX=FALSE. 
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4.  The  Limit  of  a  Quotient  of  RP-forms 

An  RP-form  Is  defined  as  an  expression  obtainable  by 
combining  polynomials  and  positive  rational  number  powers  of 
polynomials  using  the  operators  ♦,  -,  *.  Here  are  some 
examples  of  RP-forms, 

3  1/3 

f(X)  =SQRT( 4  X  ♦  5  X)  +  (X  ♦  1)  ♦  X  ♦  C 

2/3 

g(X)  *  AX  SQRTCX+l)  ♦  X 

Note  that  the  definition  of  RP-forms  does  not  allow 
nested  radicals  of  X.  For  Instance  the  following  Is  not  an 
RP-form. 


2/3 

(X  +  SQRT(X) ) 

Now  consider  the  limit  problem, 

LIMIT(E(X),X,L>,  E(X)  «N(X)/D(X) 

where  N(X)  and  D(X)/0  are  RP-forms  and  at  least  one  of  them 
Is  Irrational.  Note  that  N(X)  or  D(X)  can  be  a  constant. 

Let  us  define  the  operators  EXPO  and  COEF  by 
EXP0(N(X) ) 

“highest  exponent  of  X  In  N(X), 

C0EF( E,N(X) ) 

“the  coefficient  of  E  In  N(X). 


29 


For  example, 

EXP0(g(X) )  *7/6, 

COEF(X**(3/2),f (X))  »2 

EXPO  Is  obtained  by  a  simple  special  purpose  routine;  COEF 
Is  computed  by  using  the  RATCOEF  routine  In  the  rational 
function  package  of  MACSYMA  [16 J. 

Now  the  limit  Is  computed  by  the  following  algorithm. 

(1)  If  L  -MINF,  X  Is  replaced  by  -X,  L  by  INF. 

(2)  If  L  aINF,  the  following  asymptotic  analysis  Is 
made. 


I)  compute  a  =EXP0(N(X)),  b»EXPO(D(X ) ) . 

II)  Let  N1(X)  and  D1(X)  be  N(X)  and  D(X)  with 
polynomials  under  radicals  replaced  by  their 
leading  terms,  respectively.  Compute 

a’*EXP0(Nl),  b'=EXP0(Dl>. 


Ill)  If  either  a/a'  or  b/b',  then  It  Is  usually 
sufficient  to  replace  each  radical  In  E(X)  by 


the  first  few  terms  In  the  series  expansion  of 
the  radical  about  XaINF,  obtaining  a  rational 
function  R(X).  The  answer  Is 
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LIMIT(R(X),X,INF). 

Otherwise,  compute 

c=C0EF (X**a,Nl)/C0EF(X**b,N2).  If  b>a,  return 
the  answer  c*ZER0+  If  a>b,  return  the  answer 
c*INF,  otherwise,  return  the  answer  c. 

(3)  If  L  =  INFINITY  then 

Carry  out  steps  (2-1)  and  (2-11),  If  aj*a',  bj4b • 
or  a=b,  the  answer  Is  UNO.  Otherwise,  If  a  >  b 
the  answer  Is  INFINITY,  while  If  a  <  b  the 
answer  Is  0. 

(4)  If  L  Is  finite,  compute  N(L)  and  D(L),  then 

If  N( L)  =D(L)»0,  apply  (.'Hospital's  rule.  Since 
L'HospItal's  rule  can  not  succeed  In  case  each 
term  In  N ( X )  and  D(X)  has  a  branch  point  at  X=l., 
that  Is,  N(X)  and  D(X)  have  a  common  factor  (X- 
L)**(a/b)  with  a  and  b  positive  Integers.  In 
this  case  the  factor  Is  removed  from  N(X)  and 
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D(X)  before  using  L'Hospltal's  Rule. 

If  N(  L)  =0  and  DOJ^O, 

the  answer  Is  0  If  CPLX=TRUE,  otherwise  the  be¬ 
havior  of  D(L)*N(X)  at  L  Is  examined  to  decide 
whether  0,  ZERO*  or  ZERO-  Is  the  answer. 

If  D( L) =0  and  N(L)*0, 

the  answer  Is  INFINITY  If  CPLX  “TRUE,  otherwise 
the  behavior  of  N(L)*D(X)  at  L  Is  examined  to 
decide  whether  INF  or  MINF  Is  the  answer. 
Otherwise,  the  answer  Is  E(L). 

For  large  X  the  series  expansion 


■*■»?  *»  ^no^*'*Mi«>wv*W!erwr*yw^*^^  wt*~w**:r?rrw*i 
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5.  Indeterminate  Forms  and  L'HospItal's  Rule 

L'HospItal's  rule  Is  applied  to  an  Indeterminate  form 
E(X)  of  the  form  0/0  or  INFINITY/ INFINITY.  It  first 
calculates  N'(X)  and  D'(X),  then. 

If  D'(X)*0,  L'HospItal's  rule  falls,  otherwise,  simplify 
N'(X)/D'(X)  and  evaluate 

LIMIT  (N* (X)/D* (X),X,L) • 

The  expression  N'(X)/D'(X)  can  be  more  complex  than 
E(X)  and  successive  application  of  this  rule  may  lead 
nowhere.  Therefore  the  number  of  times  this  rule  Is  called 
successively  Is  counted  and  the  relative  complexity  of 
N'(X)/0'(X)  to  E(X)  Is  tested  to  decide  whether  to  continue 
this  approach.  0u**  criterion  of  complexity  Is  based  on  the 
number  of  distinct  nonratlonal  components  of  an  expression. 
If  this  number  grows  for  three  consecutive  times,  the  ap¬ 
plication  of  L'HospItal's  rule  Is  halted. 

The  Indeterminate  form  0*1 NF INI TY  Is  transformed  to 
either  0/(1/ INFINITY)  or  INFINITY/C1/0),  depending  on  which 
Is  simpler,  before  applying  L'HospItal’s  rule.  Other  In¬ 
determinate  forms  such  as  1** INFINITY,  INF**0  and  0**0  are 
handled  by  the  logarithmic  reduction  : 


^rrvv«>^»- 
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F(X) **G(X)  «  %E**(G(X) *LOG( F(X) ) ) . 

Theoretically  the  Indeterminate  form  (INF  -  INF)  can 
always  be  reduced  to  0/0  by  rewriting  It  as 

(I)  (1/INF  -  1/INF)/(1/!NF*INF), 

but  this  method  often  makes  the  expression  much  more 
complicated.  It  Is  useful,  though,  for  expressions 
Involving  trigonometric  functions  as  can  be  seen  In 

1  1 

L I M I T ( . ,X,0,  PLUS)  *  0. 

X  SIN(X) 

DELIMITER  transforms  the  given  expression  In  this  example  to 

S I N(X)-X 
X  SIN(X), 

then  applies  L'Hospltal’s  Rule  to  obtain  the  answer  0. 
Therefore,  the  method  (I)  Is  used  for  expressions  Involving 
trigonometric  functions.  Other  types  of  expressions  can  be 
dealt  with  more  readily  by  comparing  degrees  of  Infinity  of 
the  subexpressions  or  by  series  expansion. 
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6.  On  the  Sign  of  Infinity 

If  a  function  f ( X )  Is  discontinuous  at  X*»a,  another 
function  g(f(X))  Is  not  necessarily  discontinuous  at  a. 

For  example 

l/X  A 

LIMITUA  X  +  1)  ,X,0)«  *E 

Therefore/  discontinuity  may  be  encountered  In  the  course 
of  computing  a  limit  of  a  continuous  function.  Consider 
1/E(X)  with  E(a)°0  and  the  limit  problem 

UMIT(l/ECX),X,a), 

If  CPLX  =  TRUE/  the  answer  Is  INFINITY.  Otherwise/  a  rou¬ 
tine  named  BEHAVIOR  Is  used  to  analyze  the  behavior  of  E(X) 
near  X=a.  If  E(X)  approaches  0  from  above  or  below  as  X  ap¬ 
proaches  a,  then  the  answer  Is  INF  or  MINF/  respectively. 
Otherwise/  the  answer  Is  INFINITY* 

The  BEHAVIOR  routine  uses  differentiation  to  analyze 
the  behavior  of  a  function  near  a  point  by  Investigating  Its 
slope  or  concavity  at  the  point.  This  routine  Is  also  used 
by  other  programs.  One  of  these  Is  the  program  which  com¬ 
putes  limits  of  the  trigonometric  function  tangent.  In  case 
the  argument  of  TAN  approaches  %P\/2,  say#  It  Is  Important 
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to  know  the  direction  of  approach. 

7.  Heuristics  for  Comparing  Orders  of  Infinity 

Let  f(X)  and  g(X)  be  two  real-valued  functions  which 
become  positively  Infinite  as  X  approaches  IMF.  We  define 
the  symbol  >>  by 

Definition  :  g(X)  >>  f(X)  If 

LIMIT(f(X)/g(X),X,INF>  »0. 

Using  this  symbol  we  can  make  the  following  brief 
table. 

...  %E**%E**X>>%E*X>>X**n» log(X)>>1og( log(X) )  ... 

Many  limit  problems  with  the  variable  approaching  INF 
can  be  solved  very  efficiently  by  using  this  concept.  Some 
of  these  problems  can  be  difficult  to  solve  by  other  me¬ 
thods.  For  example,  L'HospItal  rule  falls  to  compute  the 
1 1ml  t 


LIMITC (%E+1)**X**2/^E**X,X, INF), 

while  the  answer  Is  obviously  INF.  We  Incorporate  this  con¬ 
cept  by  a  routine,  STRENGTH.  It  can  classify  the  order  of 
Infinity  of  the  argument  according  to  the  following  rules. 


^  w*wwf^?'^^  t'rtrxwmm'mr*.  rww  ;mww»mw  ~ 
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STRENGTH(c)*(c> 

STREMGTH( LOG( R(X) ) )*  LOG 
STRENGTH(A**B(X))*  EXP 
STRENGTH ( P(X ) **(m/n) )*  m*d/n 

where  c  Is  Independent  of  X,  m  and  n  are  positive  Integers, 
d  the  degree  of  the  polynomial  P(X).  Note  the  STRENGTH  of  a 
constant  Is  denoted  by  that  constant  Inside  parenthesis  and 
the  STRENGTH  of  a  polynomial  Is  a  constant.  B(X)  and  R(X) 
can  be  any  functions  In  X,  except  those  which  can  cause  the 
relevant  argument  of  STRENGTH  to  be  simplified  Into  one  of 
the  other  three  cases.  We  can  assume  A  to  be  5E  for  If  A  Is 
any  other  expression  It  Is  always  possible  to  change  the 
base  to  $E.  STRENGTH  of  a  sum  Is  the  maximum  of  the 
strengths  of  the  terms  In  the  sum. 

There  Is  a  basic  comparison  routine  which  knows  the  re¬ 
lations  between  the  strengths  of  functions, 

EXP>>a»L0G>>(b)  and  c»d  If  c>d 
where  a,  b,  c,  and  d  are  constants. 

Let  N(X)  and  OCX)  be  two  products  satisfying 


L I M I T ( N ( X ) , X , IMF)* INF 
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and 

LIMIT(D(X),X,INF)-INF 

and  no  factor  In  N(X)  or  D(X)  Is  a  sum  or  a  trigonometric 
function.  To  compare  the  relative  order  of  Infinity  of 
N(X)  and  D(X),  the  following  method  Is  used.  The  answer  1 
Indicates  N(X)>>D(X),  -1  Indicates  D(X)>>N(X),  0  no 
decision. 

Algorithm  COMPARE(N(X),D(X) )  : 

1)  Remove  from  N(X)  and  D(X)  any  common  factors. 

2)  Apply  the  STRENGTH  routine  to  each  factor  of  N(X) 
and  obtain  the  maximum  strength  SN. 

3)  Do  the  same  to  D(X)  and  obtain  Its  maximum  strength 
SD. 

4)  If  SN>>SD  the  answer  Is  1,  If  $D<<SN  the  answer  Is  - 

1. 

Otherwl se, 
l )  Let 

Nl(X)=product  of  all  factors  with  strength  SN  In 
N  ( X ) . 

Dl(X)=product  of  all  factors  with  strength  SD  In 


D(X) . 

it)  If  SN“SD»some  constant#  set 
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SN “STRENGTH (Nl)  and  SD=STREMGTH(D1) . 

If  SN>>SD#  the  answer  Is  1#  If  SN<<SD  the  answer  Is 
-1.  If  both  M1(X)  and  D1(X)  are  polynomials  In  X# 
the  answer  Is  the  leading  coefficient  of  the 
polynomial  N1(X)  -  Dl(X).  Otherwise  If  che  variable 
INDICATOR  has  the  value  1#  this  means  each  of  N(X) 
and  D(X)  Is  an  exponent  of  an  exponent  of  an  expo¬ 
nential  expression  and  they  come  as  a  result  of 
recursive  calls  to  COMPARE.  For  example  N(X)  may 
come  from  %E**$E**N(X) .  In  this  case#  the  algorithm 
goes  as  follows.  First  INDICATOR  Is  set  to  0  then 
E=UMIT(N(X)-D(X)#X# INF)  Is  computed  and  the  answer 
Is  1  If  E>0#  -1  If  E<0#  0  otherwise. 

If  SN=SD=EXP# 

The  exponents  A(X)  and  B(X)  of  M1(X)  and  D1(X) 
are  computed.  That  Is  N1(X)  =%E**A(X)  and  D1(X) 
“%E**B(X).  If  both  A(X)  and  B(X)  are  polynomials  In 
X#  the  answer  Is  the  leading  coefficient  of  the 
polynomial  A(X)  -  B(X)  If  this  polynomial  Is  not  a 
constant#  0  If  It  Is.  If  A  or  B  Is  a  sum#  A  Is  set 
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to  LIMIT(A(X)-B(X),X,INF)  and  the  answer  Is  1  If 
A=INF,  -1  If  A=MINF,  0  otherwise.  In  case  neither  A 
nor  B  Is  a  sum,  If  both  A  and  B  are  of  strength  EXP, 
set  the  variable  INDICATOR  to  1.  The  answer  Is  ob¬ 
tained  by  evaluating  C0MPARE(A1(X),B1(X)),  where 
A(X)  =  %E**Al(x)  and  B(x>  =  %E**B1(X). 

If  SN=SD»LOG, 

If  N1(X)  Is  LOG( f (X) )  and  D1(X)  Is  L0G(g(X))  then 
set  SN  to  STRENGTHCf (X) ),  SD  to  STRENGTH(g(X) ) .  Mow 
the  answer  Is  0  If  both  SN  and  SD  are  constants,  1 
If  SN>>SD,  -1  If  SN<<SD.  The  answer  Is  0  otherwise. 

We  shall  next  discuss  the  Indeterminate  form  (INF-INF) 
and  see  how  COMPARE  can  be  used  In  such  situations.  Let 

n. 

F(X)  =  £  ft  (X) 

1*1 

where 

LIMITCf.  (X ) ,X,  INF)  «INF  or  MINF 

for  1  =  1,2, ...n. 

LIMIT(F(X),X, INF)  becomes  Indeterminate  when  there  exists  I 
and  j  such  that 

L I  Ml T( f j  (X),X, INF)  =  INF 

UMlTCfj  (X),X,  INF)  =  MINF 
J 


and 


-jt  .  .p’-ny  .^V'.  .• ' r  '~  '  "'  rrv''^X'<  ■'u^r‘M?'??.  V*  ;’--'TP'-i'^u'/ 
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The  method  of  orders  of  Infinity  can  sometimes  be  ap¬ 
plied  to  this  problem.  The  algorithm  used  Is  as  follows. 

The  program  forms  two  lists  LI  and  L2  of  the  f j  *  s  which 
go  to  INF  and  MINF  respectively.  Then  from  each  of  LI  and 
L2  a  member  of  maximum  order  of  Infinity  can  be  obtained. 

Let  these  be  r(X)  In  LI  and  s(X>  In  L2.  Now  If  r(X)  »s(X) 
the  answer  Is  INF,  If  s(X)>>r(X)  the  answer  Is  MINF.  Other¬ 
wise,  the  problem  can  be  very  difficult.  However  the  pro¬ 
gram  tries  to  evaluate 

INF*LIMIT(FCX)/s(X),X,  INF), 
which  sometimes  produces  an  answer. 

Now  we  will  follow  the  major  steps  of  the  solution  of  a 
problem  In  a  more  detailed  manner.  Let 

2 

X*SQRT( X  ♦  1) 

A ( X )  *  %E 

2 

X 

B ( X )  =  %E 

Consider  the  problem  L I  Ml T( A(X) -B(X) ,X, INF) ,  First  the 
Indeterminate  form  (INF-INF)  Is  encountered  through  the  fol¬ 
lowing  steps. 


LIMIT(A(X)-B(X),X, INF) 


■'  ■  ***** i:;  W* *fcw»j-  ^m^/r'v'-  ‘ 


41 

LIMIT(A(X),X,IMF)-LIMIT(B(X),X,INF) 
LIMIT(A(X),X,INF)  =%E**INF  » I NF 
LIMIT(B(X),X/  INF)  =%E**INF  =  I  NF 
INF-INF 

Then  COMPARE  Is  applied  to  A(X)  and  B(X),  only  to 
arrive  at  no  conclusion  about  the  relative  order  of  Infinity 
of  A(X)  and  B(X). 

COMPARE  (A(X  )  ,B(X) ) 

$TRENGTH(A(X))«EXP 

?  i  RENGTHC  B(X)  )**EXP 

2  2 

COMPAREC X*SQRT(X  ♦  1),X  ) 

2 

STRENGTH(X*SQRT(X  +  1))  «  2 
2 

STRENGTH(X  )  «  2 

Now  the  problem  Is  converted  to  the  following  form  and  the 
answer  Is  INF. 

INF*(LIMIT(A(X)/B(X),X, INF)-1) 

2  2 

LIMI T(X*SQRT(X  ♦  1)  -  X  ,X,INF)  =  1/2 
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INF*UE**(l/2)-l) 

INF 


8.  Examples 

A  number  of  limit  problems  solved  by  DELIMITER  are 
Included  In  this  section.  They  are  presented  In  the  form  of 
actual  Inputs  and  outputs  of  the  MACSYMA  system.  Lines 
labelled  (Cl)  are  tnput  or  command  lines  and  (Dt)  output 
1 1 nes  or  answers . 


(Cl)  LIMIT(X**LOG(  1/X)/.  X,  INF)Q 
(Dl)  0 

(C2)  (C0S(X)-l)/(  %E**X**2  -1)Q 


(D2) 

(C3)  LIMIT  (D2/X/0)Q 
(03) 

(C4)  ( 1+A*X)**( 1/X)@ 
(D4) 

(C5)  LIMIT(D4,X,0)@ 
(D5) 


COS(X)  -  1 


2 

X 

%E  -  1 
1 
2 

1/X 

(A  X  ♦  1) 


A 

*E 


(C6)  X**2*(4*X**4+5)**(l/2)-2*X**4@ 
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(C6)  X**2*(4*X**4+5)**(l/2)-2 

2 

(D6)  X  SQRT(  4 

CC7)  LIMI T( D6,X, INF)@ 

(D7) 

(C8)  1/X-l/SI N<  X>@ 

(D8) 

(C9)  LIMIT(D8,X,0,  PLUS)@ 

(D9) 

(CIO)  2E**(X*(X**2+l)**(l/2))-2E 

2 

X  SQRT(X 

(DIO)  %E 

(Cll)  LIMIT(D1,X, INF)@ 

(Dll) 

(C12)  UE**X+X*LOG(X))/(SE**(X** 


**4@ 

4  4 

♦  5)  -  2  X 

5 
4 

1  1 
X  SIN(X) 

0 

*X**2@ 

2 

+  1)  X 
-  %E 

INF 

+l)**(l/2)+LOG(X**4+X+l) )@ 


(D12) 


X 

%E  +  X  LOG(X) 


3 

4  SQRT(X 

LOG( X  ♦  X  ♦  1)  +  %l 

(C13)  L I  Ml T(D4/  X, I NF)@ 

(D13)  0 

(C14)  l/(X**3-6*X**2+ll*X-6)@ 


+  1) 


(D14) 


1 


(CIS) 
( D 1 S ) 
(C16) 

(D16) 

(C17) 

(D17) 

(C18) 


3  2 

X  -GX  +  11  X  -  6 
LIMIT(D14,X,2,MINUS)@ 

INF 

( X  *  S  QRT (X*5)+1)/(SQRT(4*X**3+1)+X)@ 

X  SQRT(X  +  5)  ♦  1 
3 

SQRT(4  X  ♦  1)  ♦  X 
LIMIT(01G,X, INF)Q 

1 

2 

TANCX )/LOG( COS(X)  )0 


(D18) 

(C19) 


TAH(X) 

LOG(COS(X)> 

LIMIT(D18,X,$PI/2,MlNUS)@ 


(019) 


MINF 
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CHAPTER  I  I  I 

KANDERER-AN  INTRODUCTION  AND  OUTLINE 

1.  Introduction 

The  present  chapter  and  the  three  subsequent  ones  are 
devoted  to  a  complete  description  of  WANDERER/  a  heuristic 
program  for  evaluating  definite  Integrals  over  a  real  range. 

WANDERER  cannot  solve  all  definite  Integration 
problems/  but  It  can  certainly  obtain  solutions  to  a  large 
number  of  Integrals.  Many  examples  presented  are  actual 
problems  taken  out  of  graduate  text  books  [1/  5].  The  basic 
approach/  In  WANDERER/  for  evaluating  an  Integral  Is  by  com¬ 
plex  contour  Integration  and  residue  techniques.  This  Is  a 
very  general  method  particularly  useful  In  evaluation  of 
Infinite  Integrals.  Some  other  methods  available  to 
WANDERER  are:  substitution/  Integration  by  parts,  diffe¬ 
rentiation  with  respect  to  a  parameter,  pattern  recognition 
and  table  look-up,  f Ini te-to-lnf Ini te  conversion.  Introduc¬ 
tion  of  a  parameter  and  partition  and  transformation  of  the 
range  of  Integration. 
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The  user's  command  tn  MACSYMA  for  definite  Integration 
Is  the  key  word  DEFINT.  DEFINT  takes  four  arguments  as  In 

DEFINTCexp,  var,  a,  b) 

where  exp  Is  the  Integrand,  var  the  variable  of  Integration, 
a  and  b  the  lower  and  upper  limits  of  Integration.  There  Is 
no  restriction  as  to  what  types  of  Integrands  are  allowed  as 
Input,  as  long  as  they  remain  elementary  functions. 

However,  WANDERER  requires  the  Integrand  to  be  finite  In  the 
Integration  range  except  possibly  at  the  end  points  a  and  b. 

A  few  special  symbols  are  used  In  WANDERER  for  some 
branches  of  the  logarithm  function.  They  are  listed  In 
appendix  D. 

2.  Outline  of  Approach 

Before  going  Into  the  details  of  methods  and  algorithms 
for  evaluation  of  the  many  types  of  definite  Integrals  which 
will  be  discussed  In  the  next  three  chapters,  a  qutck  look 
at  the  whole  picture  with  emphasis  on  the  flow  of  control  Is 
In  order. 

WANDERER  Is  a  heuristic  program  which  computes  definite 
Integrals  by  trying  to  apply  one  or  more  of  the  methods  or 


■Sl.  Vf  y*  -s  •*•  WpY>tfl..^77  •*•**  V  «*)(?•*■*’ '  '  * <"  ■*- 
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algorithms  built  Into  It.  The  clues  as  to  which  of  the 
methods  to  apply  are  obtained  by  examining  the  range  of 
Integration  and  the  form  of  the  Integrand. 

Pig.  1  on  the  next  page  serves  as  a  simple  outline  of 
the  flow  of  control  In  WANDERER. 


,mr«y  vf  ~v  «*A 


S=success 
F»fa  l  lure 


WANDERER 


simple 


ntegrals-**  return  answer 
F  * 


normal  I zatlon 
of  Integrand  and 
Integration  range 


methods  for 
Infinite  integrals, 
applied  according  to 
range  of  Integration 

3 


methods  for 
finite  Integrals, 
Including  proper  and 
Improper  Integrals 

$ 


return  answer 


more  general  methods 
transformation,  — 
change  of  variable 


$ 


^return  answer 


r 

divergence  test 


Pig.  1 

As  Indicated  In  fig.  1.  WANDERER  can  be  said  to  have 
five  stages.  In  the  first  stage  simple  Integrals  are  com¬ 
puted  which  Include  cases  such  as  a  constant  Integrand,  an 
Integral  with  equal  upper  and  lower  limits  of  Integration 
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and.  If  the  variable  of  Integration  Is  X,  Integral  with 
Integrands  which  are  polynomials  In  X  or  %E**X,  The 
normalization  stage  will  make  sure,  by  transformation  If 
necessary,  that  (1)  the  lower  limit  Is  less  than  the  upper 
limit,  and  (2)  constant  factors.  If  any.  In  the  Integrand 
are  removed  to  be  multiplied  Into  the  final  answer  of  the 
I ntegral . 

The  third  stage  comprises  most  of  the  methods  that  wtll 
be  described  In  chapters  4  and  5  and  Is  the  work-horse  of 
the  whole  program.  Many  methods  In  this  stage  are  grouped 
according  to  the  Integration  range  In  which  they  are 
appropriate.  Finite  Integrals  of  rational  functions  are 
transformed  Into  Infinite  Integrals  which  will  be  evaluated. 
For  other  finite  Integrals,  an  effort  Is  first  made  to 
determine  whether  they  are  Improper  Integrals,  The  absolute 
divergence  of  a  finite  Improper  Integral  Is  tested  before 
any  attempt  at  evaluation.  If  the  third  stage  falls  to 
obtain  any  results,  the  given  problem  enters  the  fourth 
stage  which  contains  more  general  methods  appropriate  to 
various  types  of  Integrands.  A  transformation  or  change  of 
variable  Is  often  done  In  this  stage.  When  WANDERER  runs 
out  of  methods,  the  convergence  of  the  given  Integral  will 
be  challenged.  In  case  It  Is  divergent,  WANDERER  will  so 
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Indicate  In  the  output.  Otherwise,  Indication  of  failure 
will  be  returned.  It  may  be  argued  that  tests  for 
convergence  should  be  conducted  before  evaluation,  but  this 
Is  not  necessary  because  WANDERER  Is  so  designed  that 
whenever  tt  produces  an  answer,  the  given  Integral  Is 
convergent.  That  Is  to  say  each  Individual  method  has  Its 
own  convergence  and  divergence  conditions  built  In. 

The  limit  program  described  In  the  previous  chapter  Is 
used  whenever  a  limit  computation  Is  needed  In  WANDERER.  It 
may  often  be  required  In  changing  the  variable  of  Integra¬ 
tion,  computing  residues  and  testing  for  convergence  or 
divergence. 

3.  Outline  of  the  Computation  of  Residues 

One  of  the  most  powerful  methods  In  evaluation  of 
definite  Integrals  Is  the  use  of  contour  Integration  through 
residue  calculations.  The  Importance  of  contour  Integration 
and  the  residue  theory  In  the  sequel  warrants  a  brief 
summary  of  relevant  facts  from  complex  analysis  [5], 

(I)  The  Cauchy  Integral  Theorem:  Let  D  be  a  simply 
connected  domain  and  let  f(Z)  be  analytic  In  D.  Let  C  be  a 
closed  contour  In  D.  Then 
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I  f(Z)  dZ  -  0 
->C 

A  simply  connected  domain  Is  Intuitively  an  open  set  of 
points  with  no  holes  or  cuts  In  Its  Interior.  For  example 
the  disk  | Z |  <  0  Is  simply  connected  while  the  annulus  1  < 

| Z |  <  2  Is  not.  The  Cauchy  Integral  theorem  Is  one  of  the 
most  Important  tools  In  complex  analysis.  It  Is  by  use  of 
this  theorem  that  path  of  Integration  can  be  deformed.  More 
precisely.  If  Cl  and  C2  are  two  different  curves  In  D  lead¬ 
ing  from  the  same  starting  point  v  to  the  same  end  point  w, 
as  shown  In  fig.  2, 


and  f(Z)  Is  analytic  In  D,  then 


f  f(Z)  dZ  «  [  f(Z)  dZ. 

W  JC2 


Therefore,  the  contour  Cl  can  be  deformed  Into  C2.  The 
above  Is  an  Immediate  consequence  of  Cauchy's  Integral 
theorem  which  gives 
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L 


f(Z)  dZ 

Ci 


f(Z)  dZ 


(ID  Definition  of  Residue:  Let  f(Z)  be  analytic  In  the 
punctured  neighborhood  of  a  point  p.  Let  Cp  be  a  small 
positively  oriented  circle  with  center  at  p.  Then  the 
residue  of  f  at  p  Is  defined  as 


1 

2fll 


f(Z)  dZ 
P 


The  value  of  this  Integral  Is  the  coefficient  of  the 

term 


(Z 


-1 


p) 


In  the  Laurent  expansion  of  f  about  the  point  p,  so  this 
coefficient  can  also  be  regarded  as  the  residue  of  f  at  p 


(III)  The  Residue  Theorem:  Let  C  be  a  simple  closed 
positively  oriented  contour.  Let  D  be  a  simply  connected 
domain  containing  C  and  Its  Interior.  Let  Al,  A2,  .  .  ., 
be  points  Inside  C.  Except  for  Isolated  singularities  at 
Al/  A2/  .  .  ./  Art/  let  f(Z)  be  an  analytic  function  In  D. 
Then 


An 
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f  n 

/  f(Z)  dZ  3  2%\  £  (Residue  of  f  at  Aj ) 

-'<2  j*i 

The  algorithm  for  residue  computation  used  In  WANDERER 
will  be  described  briefly  here.  Detailed  discussion  of  this 
algorithm  can  be  found  In  sect.  6-1. 

Let  f(Z)  be  a  function  of  a  complex  variable  Z, 
analytic  everywhere  In  a  domain  D  except  for  a  number  of 
poles  In  D.  Suppose  F(Z)  can  be  written  In  the  form 

f(Z)  3  U(Z)/V(Z) 

such  that  U(Z)  Is  analytic  In  D.  This  means. that  poles  of 
f(Z)  are  zeroes  of  V(Z).  Suppose  p  Is  a  pole  of  order  m  of 
f(Z),  the  residue  of  f(Z)  at  p  fs  computed  by  the  following 
a I gorl thm, 

RESIDUE  ALGORITHM: 

If  m  3  1/  compute  as  answer 

U(p)/V'(p) 

otherwise.  If  V  Is  a  polynomial, 

(I)  Set  V  to  the  quotient  of  V(Z)/(Z  -  p)**m 

which  Is  computed  by  long  division. 

(II)  Return  the  result  computed  from 

m-1 

1  l±\  U(Z> 

(m  -  1)1  I  dZ/  WIT  Z=p 


otherwise,  (m  >  1,  V  not  a  polynomial)  compute  and 
return  as  answer  the  limit 
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Ltm 

Z-P 


,  m-1 

1  d 

i  m 

(Z  -  p)  f(Z) 

(m  -  1)1  \dZj 

DELIMITER  Is  used  In  obtaining  such  a  limit.  Methods 
for  Infinite  Integrals  are  discussed  In  the  next  chapter. 
Those  for  finite  Integrals  are  Included  In  chapter  5. 
Important  algorithms  and  algorithms  common  to  many  methods 
are  detailed  In  chapter  6. 
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CHAPTER  IV 

SYMBOLIC  COMPUTATION 
OF  INFINITE  INTEGRALS 


0.  Introduction 

The  aim  of  this  chapter  Is  to  give  a  detailed  account 
of  the  methods  employed  by  WANDERER  for  Improper  Integrals 
with  an  Infinite  range.  Contour  Integration  and  residue 
computation  play  a  very  Important  role  In  many  of  the 
methods.  The  algorithm  for  Integrals  with  rational 
Integrands#  which  will  be  discussed  first/  Is  the  most  com¬ 
plete.  Sections  are  formed  according  to  the  function  types 
of  the  Integrand  and  ordered  roughly  In  Increasing  com¬ 
plexity.  Quite  a  few  examples  are  Included/  some  of  them 
with  references  to  books  or  Integral  tables  Indicated. 

1.  Infinite  Integral  of  a  Rational  Function 

In  this  section  methods  for  evaluating  Infinite  Inte¬ 
grals  of  a  rational  function  will  be  presented.  The  range 
of  Integration  Is  from  either  0  or  minus  Infinity  to 
Infinity.  Other  possible  Infinite  ranges  such  as  (a  INF)/ 


57 


(MINF  a)#  where  a  Is  finite/  can  be  converted  to  (0  INF)  by 
a  change  of  variable  In  the  Integral. 

First  let  us  consider 


L 


R(X)  dX/ 


R( X)  *  P(X)/Q(X) 


where  P(X)  and  Q( X )  are  polynomials  In  X  over  the  field  of 
complex  numbers. 

WANDERER  requires  deg(Q(X) )-deg(P(X) )  3  2  to  Insure  the 
convergence  of  the  given  Integral.  If  R(Z)  has  no  real 
poles/  then  The  Integral  L  can  be  computed  by  evaluating  the 
contour  Integral 


J 


R(Z)  dZ 


around  a  familiar  seml-cl rcular  contour  In  the  upper  complex 
Z-plane  (fig.  1).  One  can  easily  prove  that  L  =  J  as  T 
tends  to  Infinity. 
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Therefore 

(1)  L  »  2  3PI  2; I  (  £  Res(R)  over  poles  of  R  Inside  C) 

The  residue  subroutine  outlined  In  the  prevtous  chapter  Is 
used  to  compute  the  residues  needed  In  this  formula.  If 
R(Z)  has  poles  on  the  real  axis  then  L  Is  divergent. 

However/  If  the  real  poles  of  R  are  all  simple  then  the 
Cauchy  Principal  Value  of  the  Integral  L  exists#  and  can  be 
obtained  by  Indenting  the  contour  C  at  these  singular  points 
on  the  real  axis.  For  such  Integrands  WANDERER  will  compute 
this  principal  value#  which  Is  given  by 

(2)  (P)L  3  2  3PI  %l(  £  Res(R)  over  poles  of  R  Inside  C) 

+  %?\  %\(  Res(R)  over  simple  real  poles  of  R)# 

as  an  answer  to  the  Integral  L.  (P)L  stands  for  the 
principal  value  of  L.  Whenever  the  answer  Is  a  principal 
value#  the  message  PRINCIPAL  will  be  sent  to  the  user  first. 


Here  Is  an  example  solved  by  WANDERER 
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(Cl)  (X**2+A*X+B)/(X**4+10*X**2+9)@ 

2 

X  ♦  A  X  ♦  B 

(Dl)  . 

4  2 

X  ♦  10  X  +9 


(C2)  DEFINT(D1,X,MINF,  INF)@ 

%P\  B  ♦  3  %PI 
12 


(D2) 


The  expression  Dl  has  four  simple  poles:  X  s  %\ ,  -%l, 
3  %\  and  -3  $1,  Using  formula  (1)/  only  residues  at  X=%l 
and  X=3  %\  need  be  computed.  The  residues  are 
at  X»*l  ( B-l+A  %l ) / ( 16 * %\ ) 

at  X=3  %\  “(B-9+3  A  %l )/( 48  %\) 


These  values  are  computed  by  the  residue  algorithm  given  In 
the  previous  chapter. 


The  full  algorithm  depends  on  finding  the  poles  of  R. 
How  this  Is  done  will  be  discussed  In  detail  after  consider¬ 
ing  the  next  Integral. 


K 


oo 

R(X)  dX 


If  R(Z)  has  no  real  pole  which  Is  positive  or  zero# 
this  Integral  can  be  evaluated  by  Integrating  the  following 
contour  Integral  [30] 
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i 


f(Z)  dZ, 


f(Z)  =  PLOG(-Z)  R(Z) 


where  PLOG  denotes  the  principal  branch  of  the  LOG  function, 
l.e.,  the  Imaginary  part  of  LOG(-Z)  lies  between  %P I  and  - 
%P\ ,  The  contour  C  consists  of  circular  arcs  of  radii  a,  b, 
and  two  straight  lines  joining  their  end-points  as  shown  In 
fig.  2. 


One  can  verify  that  contributions  from  the  two  circular 
arcs  of  C  vanishes  In  the  limit.  Thus, 


r®  ro  2  %PI  $1 

J  =  I  PLOG(-X)  R(  X)  d X  ♦  |  PL0G( -X  *E  )  R(X)  dX 

-'o  Jco 


=  f  PLOG(X)  R(X)  dX  -|  PLOG(X)  R(X)  dX  -  2  %PI  %\\  R(X>  dX 
«'© 

=  -  2  %PI  %\  K 


-i: 


-09 


1 
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Hence  the  value  of  the  Integral  K  Is  given  by 

(3)  K  3  -(£  Res(f)  over  poles  of  R  Inside  .C) 

Therefore,  the  value  of  K  Is  -1  times  the  sum  of 
residues  of  L0G(-Z)R(Z)  In  the  complex  Z-plane  cut  along  the 
positive  real  axis. 

It  may  be  of  some  Interest  here  to  mention  that  a 
theorem  In  |30|  gives  a  value  for  the  Integral  K  which 
differs  from  the  correct  formula  (3)  by  a  sign. 

Now  suppose  R(Z)  has  poles  which  are  real  and  positive, 
then  the  Integral  K  Is  divergent.  But  Its  Cauchy  principal 
value  exists  If  every  such  pole  Is  of  order  1.  This  value 
Is  given  by  -1  times  the  sum  of  residues  of  L0G(-Z)R(Z)  In 
the  entire  Z-plane  punctured  at  Z°0,  as  given  In 

(*»>  (P)K  a  -(  £  Res(f)  over  poles  of  R). 

K  diverges  If  R(Z)  has  a  pole  at  Za0.  The  following 
examples  were  computed  by  WANDERER.  For  ease  of  reference, 
expression  (Dl)  Is  repeated  here 
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(Dl) 

(C3) 

(D3) 


2 

X  ♦  A  X  ♦  B 


4  2 

X  ♦  10  X  +9 

DEFINT(D1,X,0,,  INF)Q 

2  %P\  B  ♦  6  L0G(3)  A  ♦  6  %P\ 

48 


The  answer  Jn  (D3)  Is  obtained  by  locating  the  four 
simple  poles  of  (Dl)  and  applying  formula  (3). 


(C4)1/(X**2+X+1)Q 

(D4) 


1 


2 

(X  +  X  ♦  1) 


(C5)DEFINT(D6,X,0,INF)@ 

2  $PI 
3  SQRTC3) 


Now  let  us  turn  to  the  problem  of  finding  the  poles  of 
R(Z).  The  SOLVE  tl6]  program  In  MACSYMA  knows  how  to  solve 
a  number  of  types  of  equations.  This  routine  Is  used  In 
obtaining  the  locations  and  multiplicities  of  the  zeros  of 
Q(x)/  the  denominator  of  R(X).  The  SOLVE  program  does  this 
by  factoring  Q(X)  over  the  Integers  and  applying  formulas  to 
each  Irreducible  factor  of  degree  less  than  5.  For  factors 
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of  higher  degree  only  those  of  the  form  aX**n+b  will  be 
solved.  This  facility  can  be  augmented  by  Including  more 
special  cases  for  polynomials  of  higher  degree.  If  the 
Integrand  Is  not  rational#  the  problem  of  finding  poles  can 
be  much  more  difficult#  although  SOLVE  can  find  some  of 
these  also.  In  any  event#  the  method  of  residues  can  be 
successful  only  If  all  relevant  poles  can  be  located. 

Let  p  be  a  root  of  Q(X)  =  0  of  multiplicity  m  which  Is 
obtained  employing  SOLVE.  It  Is  not  necessarily  true  that  p 
Is  a  pole  of  R(Z)  of  order  m.  This  Is  because  the  rational 
function  package  [10]  In  MACSYMA  does  GCD  cancellations  only 
over  the  Integers  while  R(Z)  Is  a  rational  function  over  the 
complex  numbers.  Therefore  p  may  not  be  a  pole  of  R(Z)  or  p 
may  be  a  pole  of  R(Z)  of  order  less  than  m.  Fortunately 
WANDERER  can  pretend  that  p  Is  really  a  pole  of  order  m  and 
proceed  with  the  algorithm  for  computing  residues  which  will 
produce  a  0#  If  p  Is  not  a  pole#  and  the  correct  residue#  In 
case  p  Is  a  pole  of  order  less  than  m.  More  discussion  on 
the  computation  of  residues  can  be  found  In  Sect.  6-1. 

Thus#  as  the  reader  can  easily  verify#  (D2)  of  example 
(C2)#  In  page  gf  #  Is  valid  even  If  A=0  and  B=1  (or  A=0  and 
B»9).  As  another  example#  correct  answers  were  obtained  In 
the  following  Integrals  even  If  %\ ;  the  square  root  of  -1# 
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Is  treated  as  a  simple  pole  of  the  Integrand. 


(D6) 


X  -  %\ 


(X  -  2  %I)(X  ♦  1) 


(C7)  DEFINT(D6,X,MINF,INF) 


(D7) 


( C 8 )  DEFINT(06,X,0, INF) 


2  %P\ 
3 


(08) 


%\  L0GC2)  +  %PI 
3 


A  special  case  check  Is  provided  as  an  auxiliary  to  the 
above  algorithm. 

P-1 


f. 


X  dX 

N  M 
(AX  ♦  B) 


(B/A) 


P/N  /-  P/N  +  M  -  11 


%?  I 


M  -  1 


M 

B  M  S I  N( %PI  P/N) 


for  M,  N,  P  positive  Integers,  M  >  P/N,  P  not  divisible  by  N 
and  AB  >  0. 

For  the  case  M=l,  AB<0,  the  following  Cauchy  principal 


value  Is  used 
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A  1  -  P/N  %P\  %P\  P 

-(-  -)  (--)  COT( . ) 

B  AN  N 

This  special  case  Is  not  necessary  for  the  algorithm 
but  It  Is  very  helpful  as  far  as  program  efficiency  Is 
concerned,  especially  for  large  M  or  N. 


Sometimes  substantial  computing  time  can  be  saved  by 
application  of  differentiation  techniques  In  Integration 
problems.  Consider  the  Infinite  Integral 

dX 


(5) 


J  « 


L 


CD 


2  10 
®  (X  ♦  X  ♦  K) 


If  the  algorithm  for  rational  functions  were  used  straight¬ 
forwardly,  differentiation  would  have  to  be  carried  out  18 
times  to  obtain  the  sum  of  residues  In  case  both  poles  lie 
above  the  real  axis.  But  since 


J  « 


dX 


2 

X  +  X  ♦  K 


It  Is  only  necessary  to  differentiate  a  quadratic  expression 
once  to  obtain  the  sum  of  residues  to  evaluate  the  Integral 
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/'■ 


dX 


X  ♦  X  ♦  X 

After  this  Is  obtained,  9  more  differentiations  with  respect 
to  a  parameter  K  are  needed.  It  means  a  saving  of  at  least 
8  differentiations. 


In  fact,  whenever  the  denominator  Q(X)  of  the  Integrand 
Is  of  the  form 

v(X)**n 

with  v(X)  a  polynomial  and  n  S  2,  this  method  of  diffe¬ 
rentiation  may  be  applicable  If  the  degree  of  Q(X)  Is  large 
compared  to  that  of  the  numerator  P(X). 

We  have  Included  this  technique  In  WANDERER.  Let 

n  *»  deg(P(X))f 
s  »  deg(v(X)>, 

Q(X)  =*  v(X)**m. 

Here  Is  a  brief  description  of  the  algorithm. 

If  2  s  m*s-n,  this  algorithm  Is  not  applicable, 
otherwise  (m*s-n  >  2) 

(I)  Set  r  to  the  least  Integer  2  (n+2)/s 

(II)  If  m  >  r 

(a)  dompute  the  original  Integral  with  m 


67 


replaced  by  r  and  v(X)  replaced  by  v(X)  ♦  ZP 
where  ZP  Is  a  new  parameter  Introduced  Into  the 
problem.  Let  the  answer  thus  obtained  be 
ANS(ZP) . 

(b)  Return  the  following  as  the  answer 


m-r  | 

m-r  (r-l)l  /  d  ) 

(-1)  -  - 1 

AfIS  (ZP) 

(m-1) !  \dZPj 

( iii )  otherwise  (m=r)/  the  algorithm  Is  not  necessary 
and  thus  will  not  be  applied. 


The  parameter  ZP  used  In  this  algorithm  makes  the 
symbolic  differentiation  In  step  (ll-b)  possible  at  all 
times.  The  final  answer  Is  obtained  by  setting  ZP  to  zero. 
Such  a  parameter  will  be  used  again  later  and  will  be 
referred  to  as  the  zero  parameter. 


For  example/ 

(D9) 


1 


2  3 

(X  ♦  X  ♦  1) 


(C10)DEFINT(D9/X/0/INF)@ 

4  %P\  1 


9  SQRT(3)  2 

In  this  example/  we  have  n=*0/  m=3  and  s  =  2.  Thus  r=l 
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and  the  value  of 


dX 


2 

♦  X  ♦  1) 

Is  computed  first. 
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2.  Infinite  Integrals  Involving  Algebraic  Irrational  Func¬ 
tions 


2.1  Evaluation  by  Residues 


Let  R(X)  be  a  rational  function. 

I 

t**K*R(X)  dX, 


°  ■/*" 


If 


Consider 

-1<  K  <  1 


LIMIT(X**(K+1)  R(X)  ,  X,  0,  PLUS)  -  0, 
LIMIT(X**(K+1)  R(X ) ,  X,  INF)  -  0 , 
and  R(X)  has  no  poles  of  order  greater  than  one  on  the 
positive  real  axis,  then  D  can  be  evaluated  by  applying 
residue  theory  to  the  contour  Integral 


f(Z)  dZ, 


f(Z)  «  (-Z) **K  R(Z) 


around  the  contour  C  shown  In  fig.  2.  The  limit  conditions 
given  above  are  convergence  tests  for  D.  They  are  computed 
by  DELIMITER.  In  computing  these  Integrals,  the  Integration 
program  will  first  obtain  two  sums  of  residues  SI  and  S2  of 
the  function  f(Z)  by  executing  the  following  two  steps. 


If  R(Z)  has  pnles  off  the  positive  real  axis,  set  SI 
to  £Res(f(Z))  at  these  poles,  otherwise  set  SI  to 
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0. 

If  R(Z)  has  simple  positive  poles#  set  S2  to 
£]Res(f(Z))  at  these  poles#  otherwise  set  S2  to  0. 

Now  the  answer  to  the  Integral  D  Is  given  by 


%P\  SI 

SINUPI  (K  ♦  1)) 


%PI  COT ( % P I  (K  +  1))  S2 


For  example#  MACSYMA  produced  the  following  result. 


(Cl)  l/((l+X)*X**(l/2))@ 
(Dl) 

(C2)  DEFINT(%#X#0# INF)Q 
(D2) 


1 

SQRT(X)  (X  ♦  1) 

%PI 


[3] 


2.2  Integrals  Related  to  the  Beta  Function 
From  the  definition  of  beta  function 


BETA(K#S) 


K-l) 


(1-X)**(S-1)  dX# 


where  K  >  0  and  S  >  0#  One  can  deduce  the  following  rela¬ 
tion 
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rao 

X 


K  -  1 


i 


dX 


(CX+D) 

where  CD  >  0  and  S  >  K  >  0. 


S 

BETA(K,S  -  K) 


K 

C 


By  a  further  change  of  variable  In  this  formula,  the 
following  very  useful  relation  Is  obtained. 

K  -  1 


(a) 


r"x 

Jo  ' 


dX 


BETA(A, B) 

A  B 
C  D  r 


r  >  0 


(CX  +D) 

where  A  ■  K/r  >  0,  and  B  ■  S-A  >  0. 

This  rather  general  formula  covers  many  Interesting 
Infinite  Integrals  Involving  algebraic  Irrational  functions. 
This  Is  built  Into  the  programs  by  using  special  purpose 
pattern  recognition  routines  (see  chapter  6)  to  examine 
whether  the  Integrand  Is  of  the  particular  form  (a).  A 
simplification  routine  for  Beta  fuijctlon  was  also  needed  to 
put  the  results  In  a  simpler  and  better  looking  form. 


(D3) 


K 

X 


X  ♦  3 


(C4)  DEFINT(D3,X,O,INF)0  [12] 

IS  THE  EXPRESSION 
K  ♦  1 

POSITIVE,  NEGATIVE,  OR  ZERO 
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POSITI VE@ 

IS  THE  EXPRESSION 
-  K 

POSITIVE/  NEGATIVE/  OR  ZERO 
POSITI VE@ 

K 

(04)  3  BETA(K  ♦  1/  -  K) 


2.3  Evaluation  by  Trigonometric  Substitution 


If  R(X/Y)  Is  a  rational  functions  In  the  two  variables 
X  and  Y/  then  an  Integral  of  the  form 


U 


SQRT(X**2-A**2) )  dX 


can  be  transformed  to  an  Infinite  Integral  of  a  rational 
function  which  has  been  discussed  In  Sect.  4-1.  The  trans¬ 
formation  below  Is  easy  to  verify. 

Let 


2  2 
SQRT(X  -  A  ) 
Y . 


X  ♦  A 


U  becomes 


4  A 


2 

A  (Y  ♦  1)  2  A  Y 

R( . / . ) 

2  2 
1  -  Y  1  -  Y 


Y  dY 


2  2 

(1  -  Y  ) 
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Now  let  Z  a  Y/(l  -  Y)  In  the  above  Integral ,  It  can  be  shown 
that  U  Is  equivalent  to 


A  (2  Z  ♦  2  Z  ♦  1)  2  A  Z  (Z  ♦  1) 

(I)  4  A  /  R( . , . ) 

L  2  Z  ♦  1  2  Z  ♦  1 


Z  (Z  ♦  l)dZ 
2 

(2  Z  ♦  1) 


Although  the  above  Is  all  that  Is  needed  for  the  given 
Integral  LI,  It  Is  sometimes  more  efficient  to  transform  U 
Into  a  finite  Integral  of  a  rational  function  of  the  trigo¬ 
nometric  functions.  That  Is,  by  setting  COS(t)  =  A/X  In  U, 
one  may  deduce  that 


(ID 


*  A  SIN(t) 

R( . ,A  TAN(t))  ( . )  dt 

COS(t)  2 

COS  (t) 


doth  of  these  methods  (I)  and  (ID  have  been  Included 
In  the  programs.  Our  heuristic  rule  for  applying  these 


transformations  Is  as  follows. 


Apply  (ID  If  It  transforms  the  given  Integrand  Into 
the  form 

COS(t)**m  SIN(t)**n, 
otherwise  apply  (I). 


For  example. 
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<05) 


1 


2  2  2 
X  SQRTCX  -  A  ) 


(C6)DEFINT(D5,X,0,INF)@ 

1 

(D6) 

2 

A 


2.4  Differentiation  with  Respect  to  a  Parameter 

We  have  Included  two  formulas  In  17ANDERER  for  Infinite 
Integrals  Involving  algebraic  Irrational  functions.  They 
are 


dX 


3/2 

♦  B  X  +C) 


1 

SQRT(C)  (B/2+SQRKAC)) 


with  A  3  0,  C  >  0,  and  B  >  -SQRT(AC),  and 


U(A,B,C), 


X  dX 


3/2 

+  B  X  +C) 


1 

SQRT(A)  (B/2+SQRKAC) ) 


V(A,B,C) 


where  A  >  0,  C  20  and  B  >  •  SQRT(AC). 


They  are  Included  because  any  Integral  of  the  form 
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(ill) 


M 

X  dX 


N+3/2 

+  B  X  +C) 


can  be  computed  from  these  known  results.  Our  tool  In  doing 
this  Is  differentiation  with  respect  to  an  appropriate 
parameter.  Let  Zl,  Z2  and  Z3  be  three  zero  parameters  and 


N  SQRT( %PI )/2 

H  =  (-1)  . 

3 

GAMMA (-  ♦  N) 
2 


The  algorithm  is  as  follows. 


Upon  deciding  that  the  Integrand  Is  as  given  In  (III), 
WANDERER  has  obtained  the  values  of  M,  N,  A/  B  and  C.  The 
answer  to  the  given  Integral  Is  then  computed  using  the 
simple  procedure: 

If  M  »  0  and  N  **  0,  return  U, 

If  M  a  x  and  N  °  0,  return  V , 

If  2N+2  ^  M,  the  given  Integral  Is  divergent. 

I f  N  ^  M,  return 

(  M  1 


H 


N-M  M 
dZ3  dZ2 


U(A,B+Z2,C+Z3) 
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If  M  ■  N+l,  return 


N 

H  [—4  V(  A,B+Z2,C), 

\dZ2  | 


If  2<H+1)  >  M  >  N+l/  then 

(I)  If  M  Is  even  set  r  ■  M/2,  return  the  answer 
/  \ 


H 


N 


N-r  r 
dZ3  dZl 


U(A+Z1,B,C+Z3>> 


(II)  If  M  Is  odd,  set  r  **  (M-l)/2  and  return 

N 
d 

H  ( .  V( A+Zl, B,C+Z3>) 

N-r  r 
dZ3  dZl 

‘  /. 


A  simple  example  ts 


(D7) 


(X  +  X  ♦  1) 


15 

2 


(C8)  DEF INT(D7, X,0, I NF ) Q 

1024 
6567561 


(D8 ) 
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3.  Infinite  Integrals  Involving  Trigonometric  Functions. 

3.1  Integrals  from  Minus  Infinity  to  Infinity 

Let  R(X)  be  a  rational  function  which  has  no  real  poles 

and 

LIMIT(R(X),X, INF)  «  0. 

The  Integral 


fm  ( 21  m  X) 

L  -  I  *E  R(X)  dX,  m  >  0 

J- CD 

Is  convergent.  Its  value  can  be  obtained  by  evaluating  the 
contour  Integral 


J 


<21  m  Z) 
2E 


R(Z)  dZ 


around  the  contour  C  given  In  fig.  1. 


Let  Cr  be  a  circular  arc  with  center  at  Z  «*  0,  radius  r 
and  argument  t,  tl  <  t  <  t2.  Jordan's  Lemma  [7]  shows  that 
If  f(Z)  approaches  0  uniformly  on  Cr  as  r  approaches  INF, 


LIMIT 

T-^+co 


(21  m  Z) 


f(Z)  dZB0, 


for  tl  2  0,  t2  ^  2P I 


and  that 
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LIMIT 


i-V 


m  Z) 

f(Z)  dZ-O/ 


for  tl  ^  *PI/  t2  ^  2%PI . 


Therefore  L  *  J  as  r  tends  to  Infinity  and  the  method 
of  residues  can  be  applied.  Moreover/  for  m  <  0  one  may 
use  the  same  method  by  using  a  contour  similar  to  C  In  the 
lower  half  complex  Z-plane. 


For  real  m/  n  and  p#  let  T(X)  be  SIN(mX)/  COS(nX), 
2E**(2I  P  X)  or  a  function  Involving  sums  and/or  products  of 
these  functions.  In  complex  exponential  form  T(X)  Is  a  sum 
of  constant  multiples  of  functions  of  the  form  %E**(%l  k  X ), 
k  real. 

Therefore 

r 

/  T(X>  R(X)  dX 

can  be  Integrated  using  complex  contour  Integral.  In  doing 
such  a  problem/  terms  In  the  Integrand  are  sorted  Into  two 
parts.  One  part  requires  a  contour  In  the  upper-half  Z- 
plane/  the  other  a  contour  In  the  lower-half. 

Before  discussion  of  additional  methods/  let  us  see  a 
few  examples  computed  using  methods  discussed  so  far. 
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(Dl) 

(C2)  DEF I NT(D1, X,M!  NF,  INF)Q 
(02) 


COS(X) 

2 

X  ♦  1 

*PI 

*E 


(C3)  X*SIN(X)/(X**2+1)@ 

(03) 

(C4)  DEFINT(D3,X,MINF,  INF)0 
(04) 

(C5)X*COS(X)/(X**2+l)@ 

(05) 

(C6)  DEFINT(D5/X,MINF,INF)@ 
(06) 

(C7)  X*COS (X)/(X**2+X+1)@ 
(07) 


X  SIN(X) 
2 

X  ♦  1 

%?  I 
*E 


X  COS(X) 


2 

X  ♦  1 


0 


X  COS(X) 


[12] 


[12] 


(C8)  DEF I NT(D7,X,MI NF, INF)@ 


2 

X  ♦  X  ♦  1 
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(D8) 

1  (SQRT(3)  -  %l)/2 

(%PI/2)  (C . %\)  *E 

SRRK3) 

1  (-  SQRTC3)-  %l)/2 

♦  (SI . )  %E  ) 

SQRT(3) 

(C9)  1/(?E**(%!*X)*(X**2+1))0 

1 

(09)  . 

SI  X  2 
SE  (X  ♦  1) 

(CIO)  DEFINT  (D9,X,MINF, INF)@ 

*PI 

(DIO) 

SE 

(Cll)  SIN(X)*D9@ 

SIN(X) 

(Dll)  . 

S  I  X  2 
SE  (X  ♦  1) 

(C12)  DEFINT(D11/X,MINF,INF)P 

SPI 

(D12 ) 

SE 

(C13)  D11*(U**(-5I*X)) 

SIN(X) 

(D13 )  . 

2  %l  X  2 
SE  (X  ♦  1) 


(C14)  DEFINKDIS/X/MINF^INF)® 
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(014) 


%P\ 


2 


*E 


3.2  Integrals  from  0  to  Infinity 

Let  us  now  discuss  some  Integrals  Involving  trigonome¬ 
tric  functions  over  the  range  (0  INF).  First  consider  the 
Integral  s 


and 

where  k  a 


II-/  C0S( k  X  )  dx 


-f 

•f 


n 


12  »  SINCk  X  )  dX 


nonzero  real  constant  and  n  >  1. 


Here  again,  the  residue  theory  can  be  applied  to 
evaluate  the  Integrals.  This  time  the  shape  of  the  contour 
Is  slightly  different.  The  contour  Is  a  sector  of  a  circle 
consisting  of  a  portion  of  a  circular  arc  with  center  at  the 
origin  and  two  straight  lines  Joining  the  end  points  to  the 
origin,  (fig.  3)  The  sector  angle  depends  on  n  and  Is 
%PI/(2n). 

Without  loss  of  generality,  let  us  assume  that  k  >  0. 
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•ty-3 


By  Cauchy's  Integral  theorem. 


L 


n 


SE 


(SI  k  Z  > 


dZ  -  0. 


Since  contribution  from  the  circular  arc  vanishes  as  r 
approaches  INF,  It  can  be  shown  that 


I 


AO  n 

(SI  k  Z  ) 

SE  dZ 


SE 


SI  SPI/(2  n) 


r<X>  n 
/  -k  r 
I  SE  dr 


The  Integral  In  the  right  hand  side  of  this  equation  can  be 
expressed  In  terms  of  GAMMA  function  (see  Sect.  5).  Let  us 
assume  that 


dr  *»  G 
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Therefore/ 
relation/  we 


taking  real  and  Imaginary  parts  of  the  above 
have 


and 


11  =  C0S(*PI/2n)  6 

12  »  SIN(*PI/2n)  6. 


For  example 

7/3 

(D15)  COSO  X  ) 

(CX6)  DEFINT(D15/X/0/INF)Q 


3  3  %?\ 

3  GAMMA  ( -  )  C0S( . ) 

7  14 

(016)  . 

3/7 

7*9 


(C17)  S I N( 9*X**( 7/3))Q 

7/3 

(017)  S I N( 9  X  ) 

(C18)  DEFINT(D17/X/0/INF)@ 

3  3  %?\ 

3  GAMMA ( - )  SIN( . ) 

7  14 

(D18)  . 

3/7 

7*9 


Incidentally/  If  n  5  1#  II  and  12  do  not  converge 
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A  more  general  Integral  than  II  and  12  Is 

/”®  m  n 

J  *  /  X  EXPUI  k  X  )  dX 

Jo 

where  n  >0,  R1(m)  >  -1,  k  real  and  nonzero  and  n-R1(m)  >  1. 

Vie  shall  consider  the  case  k  <  0.  The  case  k  >  0  Is 
entirely  analogous.  Let  us  take  a  sector-shaped  contour  In 
the  fourth  quadrant  as  shown  In  fig.  4. 


By  Cauchy  Integral  theorem. 


n 

<*l  k  Z  ) 

*E  dZ  **  0 


which  Impl les  that 


J  «  EXPl-SPI  %\ <m+l)/(2n))/  r"*  EXP(k  R°)  dR 


f 
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s 

3  EX  PC  —  %P  I  %\ (m+l)/(2n))  GAMMA  ( - . ) 

n(-k)S 

where  s  3  (m+l)/n. 


Let  us  look  at  an  example. 


(D19) 


(2  *1  ♦  3) 
X 


3 

($1  X  /  2) 
*E 


(C20)  DEFINT(D19/X,0,INF)@ 

2  *1  +  4 

GAMMA( . ) 

3 

(020)  . 

(2  %\  +4)  /  3 

3  (SI  /  2) 

A  similar  result  for  k  >  0  Is  given  by 


(m+l)/n 


J  3  EX P C % P 1  *1 (m+l)/(2n))  GAMMA ( . . ). 

(m+l)/n 
k  n 


The  following  Integrals  can  be  obtained  readily  from 
the  above  relation. 


: 


co  n  m 

C0S(k  X  )  X  dX 


and 


L 


00  n  m 

SI N(k  X  )  X  dX 


For  Instance/ 
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(C21)  SIN(P*X)/(X**(1/2))Q 


(D21 ) 


S I N( P  X) 
SQRT(X) 


(C22)  DEFINT(D21,X#0,INF)Q 


(  D22 ) 


SQRTUPI ) 
SQRT( 2)  SQRT(P) 


[3] 


Many  of  the  foregoing  results  depend  on  the  verifica¬ 
tion  of  the  fact  that  contributions  from  certain  parts  of  a 
contour  vanish  after  taking  a  limit.  We  hav.e  omitted  these 
proofs.  A  typical  such  proof  Is  to  show 


I  f(Z)  dZ 

Jen 


f(Z)  dZ  »  0. 


The  above  derivation  has  relied  on  this  fact.  The  proof  Is 
In  appendix  B. 


Another  Interesting  Integral  Is 


co  K  -N 

SIN  (X)  X  dX 


where  N  and  K  are  positive  Integers#  K  N  >  2  and  (K  ♦  N) 
even.  The  Integral  U  can  be  evaluated  by  use  of  the 


recurrence  relation 
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r 


r  -s 
SIN  (X)  X  dX 


r  (r  -  1) 
(s-1) ( s-2) 


r-2 


(X) 


-(s-2) 

X  dX 


♦ 


2 

r 


(s-1) (s-2) 


r  -(s-2) 

SIN  (X)  X  dX, 


where  r  >  (s-1)  >  1.  As  one  can  see/  repeated  application 
of  this  relation  will  reduce  U  to  a  sum  of  Integrals  of  the 
same  form  as  U  but  with  N  ■  1  or  2.  For  N»l,  AUDITOR  uses 
the  following  formula 


r  P  -i  rm  p-i 

I  SIN  (X)  X  dX  »|  SIN 


(X)  dX, 


p  >  0  odd. 


The  Integral  on  the  right  hand  side  can  be  evaluated  easily. 
Integrals  of  this  type  are  considered  In  Sect.  5-2.2.  For 
N=2,  we  have 


q  -2 
IN  (X)  X 


dX 


$PI  / q — 3 / 2' 
2  1  q-1  i 


where  q  Is  an  Integer  2s  2.  For  example 


(Dl) 


S I N(R  X) 
X 


,  'V?  .  <’•  1'* 


(C2)  DEF INT(5,X, 0,  INF)Q 
(02) 

CC3)  SIN(Q*X)**2/(X**2)@ 
(03) 

(C4)  DEFINT(*,X,0,INF)@ 
(04) 


%?  I 

mm  m 

2 

2 

SIN  (Q  X) 

2 

X 

*PI  Q 
2 
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'  H-'  f *!■■■■ 
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4.  Infinite  Integral  Involving  Logarithm  Functions 


4,1  Evaluation  by  Recursion  and  Contour  Integration 


Let  R(X)  be  a  rational  function  which  Is  even  In  X. 
The  Integral 


V(N) 


i 


CO  N 

LOG  (X)  R(X )  dX, 


R(X>  =  P(X)/Q(X), 


with  deg(Q(X) )-deg(P(X) )  s  it  can  be  evaluated  by  applying 
residue  theory  In  a  recursive  manner.  Recall  that  LOG  Is  an 
abbreviation  of  PLOG  when  the  argument  Is  real  and  positive. 
Consider  the  contour  Integral 


r  n 

J(N>  »/  F(Z)  dZ,  F ( Z )  =  PLOG  (Z)  R(Z) 

k 

where  C  Is  the  Indented  contour  In  fig.  5.  As  r  approaches 


IMF  and  6  approaches  0/  contributions  from  CR  and  C* vanish 
Therefore 
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J(N) 


N  %PI  V 

R( t)  LOG  (%E  t)  dt  ♦  V(N) 


N 

R(t)  (LOG(t)  ♦  *PI  %\)  dt  ♦  V(N) 


where 


K  K  HmK 

%? I  *1  LOG  (t)  dt  +2  V(N) 


Is  a  binomial  coefficient.  Thus, 


for  N  ^  0, 

J  ( N )  lr^/Nl  K  K 
V(N)  - . \  *PI  V  V(N  -  K) 

K=t 


where  J  Is  obtained  by  residues.  Therefore  the  result 
obtained  here  Is 

V(N)  a  %PI  %l(  £  Res(F)  over  poles  of  R  Inside  C) 


K 

%\  V(N  -  K). 


Based  on  this  recurrence  relation,  V(N)  can  be  com¬ 
puted.  Such  an  algorithm  has  been  Included  In  the  programs. 
The  following  Is  an  outline  of  this  algorithm. 
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1)  Declare  each  of  AV  and  AJ  to  be  an  1  X  N  array 

2)  Set  l  to  0  then  go  to  step  4 

3)  Set  AJ ( i )  to  J(!)/2 

4)  If  I  «  N,  return  the  answer  V(N) 

5)  Set  AV ( I )  to  V(l) 

6)  Set  I  to  l+l  then  go  to  step  3 

Computation  of  V(l): 

If  I  *  0/  return  the  value  of  the  ’ntegral 


GO 

R(X)  dX, 


(storing  away  Information  about  poles  and 
corresponding  residues  of  R(Z)  for  computation  of  J), 
otherwise,  compute  from  the  formula 

1  r1-'  /l\  K  K 

V(l)  «*  A  J  C  t )  -  -  Nil  %PI  v  AV(l  -  K). 
Computation  of  J C I > / 2  5 

Using  existing  Information  of  poles  and  residues  of 
R(Z)  compute  by  residue  theory  from  the  formula 


%P\  %l(  £  Res(F)  over  poles  of  R  Inside  C). 


For  example. 


.  , ... ■* "c‘  T 


>t».i.,'¥ '  :t  •• 
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(01) 


LOG  (X) 
2 

X  ♦  1 


(C2)  DEF I NT(D1,X, 0,  INF)Q 


[1] 


(D2) 


%PI 

8 


4.2  Method  of  Differentiating  and  Introducing  Parameters 

A  very  useful  method  for  Integrands  Involving  LOG(X)  Is 
differentiation.  Consider  an  Integral  of  the  form 


A(K) 


■f 


R(X)  X  LOG(X)  dX, 


K  f  0  and  -1  <  K  <  1, 


where  R(X)  Is  rational  In  X  and  k  Is  a  parameter  which 
occurs  nowhere  else  In  the  Integrand.  If  the  Integral 


BOO 


r<o 


R(X )  X  dX 


Is  convergent  and  can  be  evaluated,  then  A(K)  Is  given  by 


AOO  »  —  BOO. 
dk 


This  method  of  differentiation  Is  valid  because  BOO  Is 
convergent  and  AOO  Is  uniformly  convergent  for  every  closed 
Interval  contained  In  the  set  of  points  (K  |  K|»0  and 
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For  example 


K 


(D3) 

( C 4 >  DEF I NT(D3 
K 

(D4)  L0G( 3)  3 


X  LOG(X) 


X  ♦  3 

/X/ 0, INF)Q 

K  D 

BETA(K+1,-  K)  ♦  3  (  —  BETAU, 

Dr 

D 

- BETA(K  ♦  1 

Dr 


-  K) 


[12] 


I rB( 1+K) 
) 

r=-K 


In  the  foregoing  discussion/  K  has  been  assumed 
symbolic.  Yet/  In  solving  actual  problems  K  may  very  well 
be  a  number.  This  difficulty  can  be  overcome  by  Introducing 
a  zero  parameter  ZP.  V7e  first  replace  X**k  by  X**(k+ZP)  In 
f(X).  This  permits  us  to  proceed  as  above  and  then  diffe¬ 
rentiate  with  respect  to  ZP.  After  differentiation  the 
resulc  Is  then  evaluated  at  ZP  *  0. 


The  property  of  the  logarithm  function 

-LOG(X)  *•  LOGU/X) 

can  sometimes  be  used  In  evaluating  Integrals  of  the  form 
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U 


CO 

f(X)  LOG(X)  dx 


In  fact/  If  a  change  of  variable  Y  «  1/X  Is  made  In  U 
and  If  the  new  Integral  happens  to  be 


P 


f (Y) L0G( Y)  dY# 


then  the  value  of  the  given  Integral  U  Is  0.  For  example 

1/3  -  1/3 

(ARCTAN(X  )  +  ARCTAN(X  ))  LOG(X) 

(05)  . . 

2 

X  ♦  1 

(06)  DEF  I  NKD5/X/0/  INF)@ 

(06)  0 


4.3  Integration  by  Parts 

Another  nice  property  of  the  logarithm  functions  Is 
that  their  derivatives  are  often  simpler  than  the  function 
themselves.  Because  of  this.  Integration  by  parts  Is 
frequently  a  suitable  method  for  Integrals  Involving  them. 
If  the  Indefinite  Integral 

U(X)  »  J f(X)  dX 
can  be  obtained  and  the  Integral 
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J 


= /  LOG(g(X))  f(X)  dX 


Is  convergent,  then 

I b  U(X)g’(X)  dX 

J  »  LOGCg(X) )  U(X)  I  -I  ( . ). 

Iaja  gCX) 

The  SIN  Integration  program  [20]  In  MACSYMA  Is  used  to 
compute  the  Indefinite  Integral  U(X).  WANDERER  uses  this 
method  to  evaluate  Integrals  of  the  form 


a 

L0GC1  ♦  X  )  dX, 


a+1  >  L  >  1. 


Here  Is  an  Integral  evaluated  by  the  method  of  Integra¬ 
tion  by  parts. 


7/2 

L0G(X  ♦  1) 

(D7)  . 

2 

X 

(C8)  DEFINT(D7/X/0/lNF)@ 


(D8) 


*PI 
5  *PI 

S  I  N( . ) 


7 
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5.  Infinite  Integrals  Involving  Exponential  Functions 
5.1  Method  of  Substitution 


Integrals  with  tntegrands  which  are  rational  functions 
of  %E**(K*X),  K  real  and  non-zero#  are  relatively  easy  to 
compute.  Without  loss  of  generality#  let  us  assume  K  >  0 
and  consider 


(1) 


P 

J-CO 


R(%E**(K*X))  dx. 


If  one  makes  a  change  of  variable 
(2)  Y  -  %E**(K*X), 

then  (1)  becomes 


(1/K) 


co 

R(Y)/Y  dY. 


This  Integral  converges  If  (1)  does.  That  Is  to  say 
R(Y)  has  Y  as  a  factor  If  (1)  converges.  For  the  Integral 


(3) 


co 

R($E**(K*X))  dX 


one  can  make  a  similar  change  of  variable 
(4)  S  ♦  1  =  %E**(K*X) 

which  gives 
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(5) 


co 

R(S+1)/(S+1)  dS. 


Integrals  of  the  form  of  (3)  and  <5)  can  be  evaluated 
by  contour  Integral  and  residue  theory  as  Indicated  In  sect. 
1.  The  methods  (2)  and  (4)  can  In  general  be  applied  to  any 
Integral  of  the  form 


b 

f(g(X))  dX 

and  the  resulting  Integral  may  be  much  simpler.  For 
Instance 


1 


(Dl)  . 

X/3  5 

*E  ( . .  7) 

X/3 

*E 

(C2)  DEFINT(D1,X,0,INF)Q 

7 

3  L0G(  — ) 

12 


(D2) 


(D3) 


5 

X/4 

SE 


X/2 

9  ♦  4 


(C4)  DEFINT(D3,X,MI NF,  INF)@ 


5.2  Use  of  Contour  Integration 


A  somewhat  more  Interesting  Integral  Is 


i 


CO 

R( %E**X)  PCX) 
oo 


dX 


where  PCX)  Is  a  polynomial  and  RCX)  a  rational  function  with 
complex  coefficients  such  that 

LIMITCRC%E**X),X,INF)  *  0 

and 

UMITCR(%E**X),X,MINF)  ■  0. 


Let  us  first  determine  a  polynomial  QCZ)  with  complex 
coefficients  such  that 

C6)  QCX)  -  QCX  +  2  %?\  %\)  -  PCX). 

QCZ)  exists  and  can  be  computed  by  the  method  of 
undetermined  coefficients.  Now  consider  a  contour  Integral 


J 


RUE**Z)  QCZ)  dZ 


taken  around  a  rectangular  contour  as  shown  In  fig.  6 


As  the  absolute  value  of  X  approaches  Infinity, 
contributions  from  the  vertical  segments  of  the  rectangular 
contour  vanish.  Hence, 


fa 

■/  RU 

J-CO 


$E**X)  Q(X>  dX 


r 

-  R(S 

J-co 


R($E**X)  Q(X+2  *PI  *1)  dX. 


It  follows  from  (6)  that  l«*J.  Mow  the  problem  of  evaluating 
the  integral  I  has  been  reduced  to  finding  the  residues  of 
R( %E**Z)Q(Z)  for  0  s=  lm(Z)  <  2  $Pl .  To  do  this,  the  poles 
of  R(Z)  are  obtained  first.  If  w  Is  such  a  pole  then 

GLOG(w) 

where  GLOG  stands  for  the  branch  of  LOG  with  Imaginary  part 
between  0  and  2*%PI,  Is  a  pole  of  the  same  order  for 
R( %E**Z)*Q(Z)  Inside  the  closed  contour.  All  such  poles  can 
be  obtained  in  this  manner.  One  example  Is 


(D5) 


X 

SINH(X)  -  %\ 


(C6)  DEFINT(D5,X,MINF,INF)@ 


100 


(D6) 


%P\ 


Note  that  the  Integrand  (D5)  has  a  pole  of  order  two  at 
Z  «  %Pl  %\/2» 

5.3  The  GAMMA  Function  and  Related  Integrals 


A  very  Important  function  closely  related  to  the 
evaluation  of  Infinite  integrals  Involving  exponentials  ts 
the  GAMMA  function  generally  defined  as 


GAMMA(Z) 


r 

«*/  %E**(-t)  t**(Z-l)  dt. 


Rl(Z)  >  0. 


Also  of  use  Is  Its  logarithmic  derivative#  the  PSI 
function 


d 

PSI (Z)  »(—  GAMMA(Z) )/GAMMA(Z), 
dZ 


A  simplification  routine  for  GAMMA  function  has  been 
written  to  make  use  of  the  many  properties  of  this  function. 

From  the  definition  of  the  GAMMA  function#  one  can 
derive  the  following  very  useful  relation 


B 

A  X 

dX  a 


C 

GAMMA(a)  *E 


a 

A  B 
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where  a  »  (D  ♦  1)/B,  R1(A)  >  0 ,  RKD)  >  -1  and  RKB) 

nonzero. 

WANDERER  has  programs  designed  to  recognize  this  form 
and  return  the  result.  Of  course  these  programs  have  to 
examine  the  signs  of  the  relevant  quantities  carefully 
before  generating  an  answer. 

For  example: 


-X 

(07)  %E 

(C8)  DEFINT(D7/X,0,  INF)@ 

SQRTC  %PI ) 
2 


(08) 


5.4  Integral  Related  to  the  Laplace  Transform 


Let  f(t)  be  a  function  of  a  real  variable  t,  then  Its 
Laplace  transform  L(f(t))  Is  defined  as 


L(f(t))  ■  F(s)  ** 


r  <-t  s> 

■  f(t>  dt. 

Jo 


Many  of  such  Integrals  can  be  evaluated  by  the  pro¬ 
grams.  Here  are  some  examples  produced  by  WANDERER. 
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(09) 


1 


S  T 

SQRT(T)  SE 

(CIO)  DEFINT(D9,T,0, INF)@ 

IS  THE  EXPRESSION 
S 

POSITIVE,  NEGATIVE,  OR  ZERO 
POSITIVES 


(DIO) 


SQRT(*PI) 

SORT(S) 


(Cll)  SIN(S*X)/(%E**X)@ 


(Dll) 


SI N(S  X) 


X 


(C12)DEFINT(D11,X,0,INF)Q 


(D12) 


S 


2 

S  +  1 


When  f(t)  Involves  trigonometrical  or  hyperbolic  func¬ 
tions,  they  are  expanded  Into  exponentials  before  L(f(t))  Is 
computed.  When  LOG(t)  Is  a  factor  of  f(t)  the  method  of 
differentiation  with  respect  to  a  parameter  can  be  applted 
If  f(t)  also  has  X**k  as  a  factor,  R1(k)  >  0.  In  this  case 
a  new  function  g(t,ZP)  Is  constructed  by  replacing  LOG(t) 
with  t**ZP  where  ZP  Is  a  zero  parameter  that  has  been 
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specially  Introduced.  The  Laplace  transform  of  g(t/ZP)  will 
be  computed  first.  Suppose  G  In 


r  (-t.  S) 

G(s,ZP)  =  L(g( t/ ZP) )  *  /  g(t/ZP)  *E  dt 

has  been  obtained.  Now/  as  mentioned  before  In  sect.  4,  all 
that  needs  to  be  don’e  to  obtain  L(f(t)>  Is  to  compute 
dG/dZP.  That  Is 

L(f(t))  ■  dG/dZP  at  ZP  -  0. 

This  Is  true  for 


f(t)  «  dG/dZP  at  ZP«0, 

and  the  fact  that  there  exists  some  sufficiently  small 
closed  Interval  containing  0  In  which  the  following  Integral 
Is  uniformly  convergent: 


/' 


(-t  s)  dG 

%E  - dt. 

d  ZP 


For  example/ 


-S  T  1/3 

(D13)  %E  X  LOG(X) 


(C14)  DEFlNT(D13/T/0/INF)@ 

1  4  4/3 

GAMMA C - )  ( LOG(S) -  PSI(-)  )  /(6  S  ) 
3  3 


(D14) 
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CHAPTER  V 

INTEGRALS  OVER  A  FINITE  RANGE 

0.  Introduction 

Many  finite  Integrals  are  proper  Integrals  whose 
Indefinite  Integrals  exist  In  closed  form  and  can  be  com¬ 
puted  rather  easily.  For  such  an  Integral/  the  evaluation 
method  is  very  straight-forward.  WANDERER  simply  obtains 
the  corresponding  Indefinite  Integrals  and  then  substitutes 
the  limits  of  Integration.  The  antiderivatives  are  computed 
by  use  of  SIN.  This  method  will  be  referred  to  as  the  anti- 
derivative  method. 

A  finite  Integral 

dX 

Is  Improper  If  f(X)  becomes  Infinite  at  some  point  c, 
b  is  c  &  a.  In  order  to  avoid  having  to  spend  computation 
time  looking  for  an  unknown  number  of  singularities  of  f 
between  a  and  b,  WANDERER  will  assume  that  Integrands  of 
finite  Integrals#  with  the  exception  of  rational  Integrands/ 


can  become  Infinite  only  at  the  end  points  of  the  range  of 
Integration.  This  Is  not  a  severe  restriction  and  does  not 
decrease  the  number  of  Integrals  It  can  handle.  Thts  ts 
true  In  the  sense  that  any  given  range  can  be  subdivided  to 
conform  to  the  above  convention.  This  convention  makes  It 
easy  to  determine  whether  a  given  finite  Integral  ts 
improper.  WANDERER  simply  checks  the  value  of  the  Integrand 
at  the  limits  of  Integration.  If  the  given  Integral  Is 
Improper,  Its  divergence  Is  tested  before  any  attempt  at 
evaluation.  WANDERER  uses  a  limit  test  for  absolute 
divergence  which  Is  discussed  In  Scet.  6-6.  If  the  given 
Integral  diverges,  WANDERER  will  so  Indicate  In  the  output. 
If  the  antiderivative  can  be  computed,  then  the  answer  ts 
sometimes  obtained  by  employing  a  limiting  process  when 
substituting  the  upper  and  lower  limits  of  Integration. 

In  thts  chapter,  attention  will  be  focused  on  the 
definite  Integrals  whose  correspond Ing  Indefinite  Integrals 
are  difficult  to  compute  or  do  not  exist. 

1.  Finite  Integrals  of  Rational  Functions. 


For  Integrals  such  as 
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where  R(X)  Is  rational  In  X,  a  and  b  are  finite,  WANDERER 
computes  U  by  transforming  It  Into  an  Infinite  Integral  by  a 
change  of  variable. 

Let  us  write 


(X  -  a)  (a  ♦  Y) 

(1)  Y  «  . . ,  X . , 

(b  -  X)  (1  ♦  Y) 


(b  -  a)  dY 

dX . . 

2 

(1  ♦  Y) 


then  It  Is  evident  that 


JrCO 

r  b  Y  ♦  a  dY 

R( . , . — >. 

„  Y  ♦  1  2 


(Y  ♦  1) 

This  Integral  can  be  Integrated  readily  by  methods  of  con¬ 
tour  Integration  and  other  means  discussed  In  sect.  4-1. 


As  an  example  let  R(X)  be  the  expression  (Dl) 

1 

(Dl)  - 

2 

X  -  3 

The  Indefinite  Integral  of  R(X)  from  0  to  1  computed  by 
the  substitution  method  above  Is  given  by 
(C2)  DEFINT(D1,X, 0,1)0 
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LOG( 2  -  SQRT( 3) ) 

(D2)  . 

2  SQRT( 3) 

One  knows  that  the  Indefinite  Integral  of  (Dl)  exists. 
In  fact  the  following  has  been  obtained  using  the  command 
INTEGRATE  In  MACSYMA 
(C3) INTEGRATE(D1,X)Q 

(D3)  ( LOG( 2  X  -  2  SQRT(3))  -  LOG(2  X  +  2  SaRT<3)>) 

/(2  SQRT( 3) ) 

The  reader  may  easily  obtain  a  result  equivalent  to 
(D2)  by  Substituting  In  (D3)  the  limits  of  Integration. 


2.  Rational  Functions  of  Trigonometric  Functions 
2.1  A  Typical  Application  of  Contour  Integration 


If  R(X,Y)  Is  a  rational  function  In  two  variables  X  and 
Y ,  an  Integral  In  the  form 


rZK 

U  =  /  R(C0S (X) ,S i N(X) )  dX 
Jo 

Is  easily  transformed  to  an  Integral  around  a  closed  con¬ 
tour.  By  setting 


COS(X) 


2 

Z  ♦  1 
2  Z 


Z  -  1 

SIN(X)  . . 

2  31  Z 
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dZ 

and  dX  **  ----# 

%\  Z 

the  given  Integral  U  becomes 


1  f  Z  ♦  1 

~  . 

tiJci  2  z 


2 

Z  -  1  dZ 

/“ - )  — 

2  %\  Z  Z 


where  Cl  Is  the  positively  orientated  unit  circle  with 
center  at  Z  *  0.  This  contour  Integral  can  then  be 
evaluated  by  finding  the  sum  of  residues  Inside  the  circle. 


Actually  this  transformation  can  be  applied  In  general 
to  any  Integral  In  the  form 


/. 


IK 

R( %E**(%I  X))  dX 


where  R  Is  rational#  by  the  change  of  variable 


Z  =  X) 


The  transformation  process  Is  simple.  The  key  point  In 
this  algorithm  Is  the  determination  of  whether  the  Integrand 
Is  In  fact  a  rational  function  of  %E**(%l  X). 


In  calculating  the  sum  of  the  residues#  only  poles 
Inside  the  contour  contribute.  Simple  poles  on  the  unit 
circle#  l.e.  those  with  absolute  value  1#  cause  the 
principal  value  to  be  computed.  The  Integral  Is  divergent 
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If  poles  of  order  greater  than  1  happen  to  be  on  the  circle. 
Examples: 

(C4)  C0S(X)**2-SIN(X)Q 

2 

(D4 )  COS  (X)  -  SIN(X) 

(C5)  DEFINT(D4,X,0,2**PI)Q 
(05)  *PI 

(C6)  SE**(2*%l*X)/UE**(2l*X)+3)@ 

2  %\  X 
*E 

(D6)  . 

%\  X 

*E  ♦  3 

(C7)  DEFINT(D6,X,O,2*2PI)0 
(07)  2  %?\ 

2.2  Utilization  of  the  Periodicity  of  the  Trl gonomef r  Ic 
Functions 

In  this  section  Integrals  of  functions  Involving  trigono¬ 
metric  functions  over  a  variety  of  ranges  will  be 
considered. 

Let  T  be  a  function  of  X  defined  by 

M  N 

T ( X )  «  COS  (X)  SIN  (X) 
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where  R 1  CM)  >  -1  and  R1(N)  >  -1.  The  following  formula  can 
be  deduced  from  the  definition  of  Beta  function.  (Sect.  4- 
2.2) 


(1) 


/. 


7t/2 

T(X)  dX 


M  ♦  1  N  ♦  1 

GAMMA  ( . )  GAMMAC . ) 

2  2 

2  GAMMAUN  +  M  +  2)/2) 


A  s Imple  example  Is 


(08) 


1/3  1/2 

COS(X)  SIN  (X) 


(C9)  DEFINT(D8,X,0,%PI/2)@ 


(D9) 


2  3 
6  GAMMA(-)  GAMMA(-) 

3  4 


5  GAMMA(  — ) 
12 


The  usefulness  of  (1)  Is  Increased  by  the  fact  that  It 
Is  possible  to  express  definite  Integrals  of  T(X)  over  a 
variety  of  ranges  In  terms  of  that  of  T(X)  over  (0  >PI/2). 
For  Instance#  the  following  relations  are  true  for  any  func • 
tlon  f. 


L 


f (S IN(X)#C0S(X) )  dX 


L 


Ti/Z 

(f (SIN(X)/C0S(X))  ♦  f(SIM(X)#-COS(X)))  dX . 


Ill 


l 


•  33/2 

f (S IN(X) / COS(X) )  dx 


L 


L 


'ttjZ 


f(SlN(X)/COS(X))  dX+  f(-SIN(X) / -COS(X) )  dX. 


More  generally/  let  f(X)  be  a  periodic  function  of  X 
with  period  2%P I .  That  Is 

f(X  +  2%P I )  •»  f(X). 

An  Integral  of  f(X)  over  some  range  (a  b)# 


S  = 


/ 

-'a 


f(X)  dX/ 


can  always  be  vr;<r^n  as  a  sum  of  Integrals  In  the  form 


(2) 


r** 

nLf‘ 


f(X>  dX  ♦  /  f(X)  dX 


f  f(X)  dX  -  r 

Jo 


f(X)  dx 


for  some  Integer  n  and  2%?\  >  c,  d  S  0.  This  Is  true  for 
there  exist  Integers  p  and  q  such  that 


a  ■  2  p  fcPI  +  d/ 
and  b  a  2  q  %PI  ♦  c. 

Then  S  Is  equivalent  to  the  sum  (2)  with  n  =  (q  -  p). 

Programs  have  been  written  to  perform  this  reduction 

and  they  are  applied  when  the  Integrand  has  a  period  2?sPI 
and  the  difference  (a  -  b)  has  %PI  as  a  factor,  some 
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examples  computed  by  WANDERER  are 

1/3  2 

(DIO)  COS(X)  SIN  (X) 

(Cll)  DEFINT(D10,X,-%PI/2,2PI/2)@ 


(Dll) 


18  SQRT(%PI)  GAMMA(-) 

3 


7  GAMMA(-) 
6 


(C 12 )  C0S(X)**3*SIN(X)**2Q 


(D12) 


3  2 

COS  (X)  SIN  (X) 


(C13)  DEFINT(D12,X,3*£PI/2,3nPI)@ 

2 

(D13) 

15 


3.  Finite  Integrals  of  Algebraic  Irrational  Functions 

3.1  Rationalizing  the  Integrand 

If  R(X,Y)  Is  a  rational  function  In  X  and  Y,  the  Inte¬ 
gral 

fC  2  2 

(1)  K  =  I  R(SQRT( A  -  X  ),X)  dX, 

Jb 

where  A  Sc  >  b  ^ -A,  can  be  rational  I  zed. 
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Let  us  write 


(2)  Y 
which  gives 

(3)  X 


2  2 
A  -  SQRT(A  -  X  ) 


X 


2  A  Y 


2 

1  ♦  Y 


v(X), 


Substituting  (3)  for  X,  (1)  becomes 


K  »  2  A 


rvcc)  Ad  -  y  ) 

J  RC . , 

-'tr<b>  i  .  v2 


2  A  Y 


1  +  Y  1  +  Y 


2 

(1-Y  )  dY 


2  2 

(1  +  Y  ) 


which  is  an  Integral  of  a  rational  function.  Integrals  of 
this  type  have  been  discussed  In  Sect.  1. 


Similarly  for  the  Integral 


/; 


2  2 

R(X,SQRT(X  -  A  ))  dX 


where  c  >  b  3A  >  0/  the  change  of  variable 

2  2 

SQRT( X  -  A  ) 

t  =  .  =  u  (X ) 

X  ♦  A 


can  be  made  to  convert  J  to 


114 


4  A 


+  1)  2  A  t 
------ 

2  2 
t  1  -  t 


t  dt 
2  2 


(1  -  t  ) 


which  Is  also  an  Integral  of  a  rational  function, 
two  such  Integrals  evaluated  by  WANDERER. 


(Dl) 

(C2)  DEFINT(D1,X,3,4)@ 
(D2) 


1 


2 

X  SQRT(X  -  9) 


%\  3  ♦  %l  SQRT(7) 

—  PLOGC - r> 

3  4 


(C3)  1/((X+1)*(4-X**2)**(1/2))Q 


(D3) 


1 


2 

(X  +  1)  SQRT(  4  -  X  ) 


Here  are 


[12] 


(C4)  TEFINT(D3/X/O/2)0 


(04) 


L0GC2  +  SQRT( 3) ) 
SQRT( 3) 


Another  method  Is  to  try  to  transform  the  given  Inte¬ 
gral  to  an  Infinite  Integral  by  the  change  of  variable  given 
In  (1)  of  Sect.  1.  For  Instance  the  Integral 

a 

R(X,SQRT((X  -  A) (B  -  X)))  dX 
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can  be  converted  to 


(B  ♦  A)  Y  (B  -  A)  Y 
. , . ) 


Y  ♦  1  Y  ♦  1 


<iY 


(Y  +  1) 


2 


The  method  of  rationalization  can  also  be  applied  to 


(CX  +  D) **( 1/Q) )  dX 


where  Q  Is  an  Integer/  C  and  D  are  constants,  and  the  range 
of  Integration  needs  not  be  finite. 


For  Integrals  of  this  type  the  substitution 

1/Q 

Y  ■  (C  X  ♦  D)  , 


Q 

Y  -  D 

X  = . / 

C 

Q  -  1 

Q  Y  dY 

dX . 

C 


will  convert  the  given  Integral  Into  that  of  a  rational 
function  In  Y. 
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The  principal  task  In  this  conversion  method  Is  to 
recognize  that  the  given  Integrand  Is  of  the  particular  form 
and  to  obtain  the  Integer  Q.  In  essence,  the  algorithm  for 
dotng  this  ts 

(I)  Obtain  a  list  L  of  all  dlsttnct  Irrational 
parts  In  the  given  Integrand. 

(II)  If  elements  In  L  are  fractional  powers  of  an 
Identical  linear  polynomial  In  X,  the  pattern  Is 
matched  and  Q  Is  set  to  the  1cm  of  the  denominators 
of  all  the  exponents  of  the  linear  polynomial, 
otherwise  the  pattern  Is  not  matched. 


3.2  Integrals  Related  to  the  BETA  Function 


The  BETA  function  Is  defined  by  the  Integral 


BETA( 


rl  k  -  i 

K,L)  =  /  X 


L  -  1 

(1  -  X  )  dX 


for  Rl(K)  >  0  and  R1(R)  >  0.  From  this  definition  one  may 
readily  deduce  the  relation 


(3) 


rv  ■  ■ 

•s  O 


C  L  -  1  1  K 

(1  -  X  )  dX  =  -  BETA(-,L). 

C  C 


WANDERER  applies  this  formula  by  recognizing  the  form 
of  the  given  Integrand.  In  chapter  6  some  techniques  and 
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programs  for  pattern  recognition  are  discussed.  For  Inte¬ 
gral  j  In  the  form 

B 

K  -  1  C  CL-1 
(X  -  A)  (B  -  X  )  dX, 

a  simply  substitution# 

Y  =  (X  -  A)/B, 
will  transform  It  Into  (3). 


118 


4.  Finite  Integrals  Involving  Logarithm  Functions 


When  the  Integrand  Is  a  function  of  LOG(X),  a  given 
Integral  may.  In  many  cases,  be  evaluated  by  transforming  It 
Into  one  which  Involves  exponential  functions.  Consider 


U 


f (LOG(X) )  dX, 


b  >  a  3  0, 


the  substitution 
(I) 


or 

(It) 


converts  U  to 


or 


I 


f(-Y)  %l 


f (Y)  %l 


-Y 

X  »  %E 


Y 

X  -  %E 


Y 


<*Y, 


dY, 


for  1  >  a, 

for  a  ^  1, 


r-  -LOG(b),  s 3  -LOG(a), 


raL0G( a),  s»  LOG(b), 


respectively.  In  case  a  =  0  and  b  =  1  (or  INF),  the  use  of 
(I)  or  (II)  will  result  In  an  Infinite  Integral  which  can 
often  be  evaluated  readily  using  methods  provided  In 
WANDERER.  For  Instance  by  use  of  (I)  the  Integral 


I. 


i  q 

L0G(X  ) 


dX 
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becomes 


q  -Y 
Y)  $E  dY, 


which  Is  an  Infinite  Integral  already  studied  (Sect.  4-5.3), 
The  method  (I)  or  (tl)  can  also  be  applied  to  Integrals  of 


the  form 


as  In 


/»' 


f (LOG(X))  dX, 


f 1  S  1 

/  LOG  (- 
Jo  X 


)  X  dX 


Jr<» 


-  R  Y  -  1 


For  Integrands  Involving  LOG(f(X)),  It  may  sometimes  be 
possible  to  simplify  the  given  Integral  by  the  substitution 

Y  =  f(X),  X  «  g(Y) 

where  the  function  g  Is  the  Inverse  function  of  f.  A  very 
simple  application  of  this  method  Is  conversion  of 


X  LOG(SQRT(X)+a)  dX 


to  the  Integral 


/: 


2  /  (Y-a)  LOG(Y)  dY. 


Techniques  of  differentiation  wtth  respect  to  a 
parameter,  as  detailed  In  Sect.  4-4.2,  can  be  employed  for 
finite  Integrals  with  a  factor 
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K 

X  LOG(X), 

In  the  Integrand. 

Examples  : 

K 

(Dl)  LOG  (X) 

(C2 )  DEFINT(D1,X,0,1)@ 

K 

(D2)  (-1)  GAMMA( K  ♦  1) 

(C3)  LOG(X)**(l/2)/(X**2)@ 

SQRT( LOG( X) ) 

(D3)  . 

2 

X 

(C4)  DEF I NT(D3,  X,  1, INF)@ 

SQRTUPI) 

(D4)  . 

2 

(C5)  X** ( 1/3)* (-LOG(X>)**( -1/2)0 

1/3 

X 

(D5 )  . 

SQRT(-l)  SQRT(LOGCX)) 

(C6)  DEF I NT(D5,X,  0, 1)@ 

SQRT( %? I  ) 

(D6)  . 

4 

SQRT(-) 

3 


RHK)  f  0 , 


[12] 


[12] 


(C7)  L0G(X)*((l-X**(l/2))/X)**(  1/2)0 


(D7) 


SQRT(1  -  SQRT(X) )  LOG(X) 
SQRT(X) 


(C8)  DEFINT(D7,X,0,1)@ 


0  3 

4  (  —  BETACr,-)) 
Dr  2 


r«*l 


(D8) 
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CHAPTER  VI 

DESCRIPTION  OF  ALGORITHMS 

1.  Computation  of  Residues 
1.1  Residue  at  a  Pole 

For  convenience  of  reference  the  algorithm  for  computa¬ 
tion  of  residues  Is  repeated  here. 

Let  f(Z)  be  a  function  of  a  complex  variable  Z# 
analytic  everywhere  In  a  domain  D  except  for  a  number  of 
poles  In  D.  Suppose  F(Z)  can  be  written  In  the  form 

f(Z)  *  U(Z)/V(Z) 

such  that  U(Z)  Is  analytic  In  D.  This  means  that  poles  of 
f(Z)  are  zeros  of  V(Z).  Suppose  p  Is  a  pole  of  order  m  of 
f(Z),  the  residue  of  f(Z)  at  p  Is  computed  by  the  following 
algorithm. 


RESIDUE  ALGORITHM: 

If  m  a  1,  compute  as  the  answer 

U(p)/V'(p) 

otherwise#  If  V  Is  a  polynomial# 
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(I)  Set  V  to  the  quotient  of  V(Z)/(Z  -  p)**m 
which  Is  computed  by  long  division. 

(II)  Return  the  result  computed  from 


1 

m-1 

d_ 

U<  Z) 

(m  -  1)1 

dZj 

V(Z) 

otherwise#  (m  >  1#  V  not  a  polynomial)  compute  and 
return  as  answer  the  limit 


m-1 

1  Id  \  m 

LIm  -  —  (z  -  p)  f(Z) 

Z-P  (m  -  1)1  I  dZl 


1.2  Evaluation  of  Contour  Integrals  by  Residue  Theory 

In  order  to  evaluate  an  Integral  of  f(Z)  around  a 
closed  contour  C  by  residue  theory#  It  Is  necessary  to 
locate  all  poles  of  f  tnslde  C.  After  this  Is  done#  the 
remaining  problem  Is  to  compute  the  sum  of  residues  of  f  at 
these  points  efficiently.  The  difficult  part  Is  finding 
poles.  WANDERER  employs  the  SOLVE  routine  In  MACSYMA  to 
solve  V(Z)=0. 

For  V(Z)  a  polynomial  In  Z#  SOLVE  finds  Its  zeros  by 


factoring  over  the  Integers  and  applying  formulas  to  each 
factor  of  degree  less  than  5.  For  factors  of  higher  degree 
only  those  of  the  form  a*Z**n+b  will  be  solved.  The  problem 
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of  factoring  large  arbitrary  polynomials  Is  non-trlvlal,  to 
say  the  least.  The  development  of  a  more  powerful  factoring 
algorithm  which  will  factor  polynomials  over  a  larger  ring 
than  the  Integers  would  certainly  be  helpful  to  SOLVE  and 
WANDERER.  The  location  of  zeros  Is  usually  more  difficult  If 
V(Z)  Is  not  a  polynomial  In  Z  but  some  more  complicated 
function.  In  such  a  case#  WANDERER  usually  uses  methods 
other  than  the  residue  theory.  An  exception  Is  when  V(Z)  is 
a  polynomial  In  £E**Z  which  has  been  discussed  In  Sect.  4- 
5.2. 

Thus  within  the  limitations  of  SOLVE/  V(Z)  will  be 
solved  and  Its  zeros  sorted  Into  a  list  of  pairs.  Each  pair 
containing  a  zero  and  Its  multiplicity/  such  as 

L  3  (  (Zl/ml)  /  (Z2/m2)  #  (Z3/m3)  /  ...). 

L  Is  then  sorted  Into  two  lists  LI  and  L2/  discarding 
poles  outside  the  closed  contour  C/  such  that 

Ll  3  a  list  of  all  simple  poles 

L2  3  a  list  of  other  poles  paired  with  their 

multlpl f cl  ties. 

At  this  point/  we  can  apply  the  RESIDUE  ALGORITHM  to 
obtain  the  desired  sum.  To  avoid  repeated  calculation  of 
V'(Z)/  the  program  checks  whether  Ll  Is  empty.  If  Ll  Is  not 


an  empty  list#  V'(Z)  will  be  computed  and  stored  for 
possibly  repeated  reference  later  In  the  computation 
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Note  that  It  has  been  assumed  from  the  beginning  of 
this  section  that  a  pole  of  f(Z)  would  be  a  zero  of  V(Z). 
This  assumption  Is  quite  reasonable  for  almost  all  of  the 
applications  In  evaluation  of  definite  Integrals  of 
elementary  functions.  However  a  zero  of  V(Z)  need  not  be  a 
pole  of  f(Z).  For  Instance  Z  3  %\  Is  a  zero  of 

4 

Z  -  1 

but  not  a  pole  of 

(Z  +  21) 

4 

(Z  -  1) 

For  each  zero  p  of  V(Z)  we  may  check  the  value  of  U(p) 
to  see  If  p  Is  really  a  pole  of  f(z).  Although  It  Is  not 
clear  what  can  be  done  If  U(p)  3  0,  since  p  may  still  be  a 
pole  of  lower  order.  A  better  method  Is  to  Ignore  the  fact 
that  If  V(p)  3  0/  U(p)  may  also  be  0  and  pretend  that  p  Is 
an  actual  pole  of  f(Z).  This  Is  valid  because  a  residue  at 
any  removable  singular  point  will  turn  out  to  be  0. 
Furthermore,  the  restdue  obtained  at  a  pole  of  order  m  Is 
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not  changed  If  any  Integer  n  >  in  Is  used  as  the  order  of  p 
In  calculating  the  residue. 

It  Is  conceivable  that  V(Z)  may  have  an  Infinite  number 
of  zeros.  Not  being  able  to  sum  Infinite  series,  WANDERER 
can  not  evaluate  Integrals  which  require  such  a  computation. 
Sometimes  only  a  finite  number  of  poles  are  tnslde  the 
closed  contour.  One  such  case  WANDERER  handles  Is  when  V(Z) 
Is  a  pol ynoml al  In  %E **Z  as  described  In  Sect.  4-5.2. 
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2.  Obtaining  Real  and  Imaginary  Parts 

In  the  course  of  evaluating  a  definite  Integral  by  com¬ 
plex  contour  Integration/  the  need  to  take  the  real  or  the 
Imaginary  part  of  an  expression  often  arise.  For  Instance/ 
to  see  whether  a  pole#  p/  lies  above  or  below  the  real  axis 
the  sign  of  lm(p)/  the  Imaginary  part  of  p#  Is  examined.  To 
determine  If  p  lies  Inside  the  unit  circle  at  Z°0/  It  Is 
needed  to  compute  ABS(p)  which  Involves  taking  the  real  and 
Imaginary  parts  of  p. 

The  algorithm  for  obtaining  Rl(p)  Is  presented  as  a 
representative  of  similar  procedures  used. 

Algorithm  REALPART(p)  ; 

1)  If  p  Is  return  0/  If  p  Is  a  number  or  any 

other  atomic  symbol/  return  p. 

2)  If  p  Is  a  sum#  (p  a£  p j )  then  return 

A 

r  REALPART(pj ) 

3)  If  p  Is  a  product  (p  a  pl*p2)  then  return 

REALPART(pl)*REALPART(  p2) 

-  IMPART(pl)* IMPART(p2! 

4)  If  p  3  %E**pl  return 

%E**REALPART(pl)*CO$( IMPART(pl)) 
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/ 


5)  If  p  3  pl**p2/  convert  p  to  $E**(p2*L0G( pi) )  then 
go  to  step  4. 

6)  If  p  =  LOG(pl),  return  L0G(ABS(pl)> 

7)  Otherwise,  return  the  form  Rl(p). 


3.  A  Heuristic  Patt’ern  Recognition  Program 

It  Is  often  the  case  that  some  pattern  recognition  Is 
needed#  at  one  stage  or  another  during  the  evaluation  of  an 
Integral.  Although  many  Integrals  can  be  evaluated  without 
any  pattern  recognition/  this  capability  remains  Important 
to  WANDERER.  One  specific  pattern  shall  be  discussed  as  a 
representative  of  such  methods  In  WANDERER.  Consider  the 
pattern 

N  M 

P(X)  ■  (B  X  +  A) 

where  B#  N/  A  and  M  are  free  of  X  (I.e.  do  not  Involve  X) 
and  all  except  A  must  be  non-zero.  This  pattern  should 
match/  for  example/  every  one  of  the  expressions 

1/2  -1/3 

X/  2  X  +  1/  (X  +  %\) 

p  q  2 

(X  -  1)  ,  and  X  +2  X  ♦  1^ 

The  last  expression  Is  the  expansion  of 
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2 

(X  ♦  1)  . 

Expressions  In  expanded  form  present  some  trouble  for 
the  recognition  algorithm.  This  difficulty  Is  overcome  by 
the  use  of  differentiation  and  rational  simplification. 

That  Is  If  an  expression  E(X)  Is  equivalent  to  an  expression 
matching  the  pattern  P(X)#  we  can  compute 

E *  ( X )  /  E ( X )  «  U(X>  /  V(X), 

cancelling  all  common  factors  In  the  numerator  U(X)  and 
denominator  V(X).  Then  V(X)  should  match  p(X)  with  M  *  1. 

As  a  result  of  this  match#  some  values  are  assigned  to  the 
variables  A  and  B.  By  use  of  these  values  the  correct  value 
of  M  and  N  can  be  recovered  from  U(X).  The  values  for  A  and 
B  thus  obtained  may  differ  from  the  true  values  by  a 
constant  factor.  This  happens  whenever  these  true  values 
have  common  factors.  The  real  values  of  A  and  B  can  be 
determined  by  comparing  A**M  to  E(0).  This  procedure  was 
suggested  by  Moses. 

The  full  algorithm  used  for  matching  PCX)  will  be 
described.  Let  It  be  called  PM.  PM  uses  another  routine  PN 
which  recognizes  the  pattern 

N 

PCX)  =  (B  X  +  A) 

As  one  can  see  In  the  following  algorithm#  the  values 
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of  B,  N,  A  and  M  will  be  set  as  the  matching  process 
proceeds.  If  the  pattern  Is  matched/  the  values  of  these 
four  variables  are  found. 

Algorithm  PM(E(X)/X)  : 

1)  If  E  8  X/  pattern  matched.  (B  =  l,N*»l,Aa0,Mal) 

2)  If  E  Is  free  of  X  or  Involves  any  of  the  func¬ 
tions:  SIN,  COS,  TAN/ LOG/ EXP,  etc.,  P(X)  Is  not 
matched . 

3)  If  E  Is  In  the  form  r**s  and  s  Is  free  of  X  then 
match  M  to  s,  otherwise  match  M  to  1.  Then,  If 
PN(r,X)  succeeds  In  matching,  the  pattern  Is 
matched. 

4)  I)  rationally  simplify  E * ( X >  /  E(X)  and  set  r 
to  the  result  obtained. 

II)  set  s  to  the  denominator  of  r,  set  r  to 
the  numerator  of  r. 

III)  If  PN(s,X)  succeeds  In  matching,  simplify 
the  expression  r  /  (N  B  X**(N-1)),  set  r  to  the 
result  thus  obtained  and  go  to  step  5. 

Iv)  P(X)  Is  not  matclv  d. 

5)  If  r  Is  not  free  of  X,  P(X)  Is  not  matched. 

6)  I)  Match  M  to  r 

II)  compute  and  set  r  to  E(0)  /  A**M 

III)  If  r  5  1,  P(X)  Is  matched.  If  r^l,  first 
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set 

1/m  1/m 

A  **  A  r  and  B  «  B  r, 

then  P(X)  Is  matched. 

Algorithm  PN<E(X),X> 

1)  If  E  Is  free  of  X,  the  pattern  not  matched. 

2)  Match  A  to  E(0)  and  set  E  to  E  -  A. 

3)  If  E  Is  of  the  form  r*X**s,  the  pattern  Is 
matched  (B  to  r ,  N  to  s).  Otherwise  there  Is  no 
match. 
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4.  A  Procedure  for  Change  of  Variables 


Substitution  of  a  new  variable  Y  for  a  subexpression/ 
say  g(X)/  of  the  Integrand  In  a  given  Integral  Is  a 
frequently  used  method  In  Integration.  Let  the  given  Inte¬ 
gral  be 


J 


f(X)  dX, 


then  the  transformed  Integral  would  be  In  the  form 


J 


r4 

»(  F(Y)  dY 


where  F(Y)  might  he  considerably  simpler  than  f(X),  F(Y)/  c 

and  d  are  computed  by  a  procedure  which  Is  called  whenever  a 
change  of  variable  Is  needed.  It  makes  use  of  two  other 
modu’  *'«s  of  MACSYMA/  namely  SOLVE  and  DELIMITER. 


I)  Use  SOLVE  to  solve  for  X  In  Y  »  g(X)/  obtaining 
X  *  h(Y)/  the  Inverse  function  of  g. 

II)  If  h  can  not  be  obtained/  return  Indication  of 
fal lure. 

til)  Compute/  using  the  DELIMITER/  c  and  d  as  In 


c  «  LIMIT(g(X),  X,  a,  PLUS), 
d  »  LIMIT(g(X),  X,  b,  MINUS) 


Iv)  Obtain  f(Y)  by  assigning  It  a  value  computed 
from 

f (h(Y) )  h'CY). 

5.  Solving  Systems  of  Linear  Algebraic  Equations 

In  sect.  4-5,2  the  need  to  compute  a  polynomial  Q(X) 
from  a  gtven  one  P(X)  satisfying  a  given  relation 

(1)  Q(X)-Q(X  ♦  2  *PI  *1)  »  P(X) 

has  been  mentioned. 

The  method  of  undetermined  coefficients  Is  used  to 
determine  Q(X),  Let  Q(X)  be  a  polynomial  In  X  with  degree 
one  higher  than  PCX)  and  unknown  coefficients  CO,  Cl,  ...  , 
Cn.  That  Is 

n+1  n 

Q(X)  *  X  ♦  Cn  X  ♦  .  .  .  ♦  Cl  X  ♦  CO. 

Equation  (1)  will  give  n  linear  relations  among  these 
coefficients.  Therefore  the  value  of  these  C's  can  be 
obtatned  by  solving  the  system  of  linear  equations  they  must 
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sat  1 sf y. 

For  the  purpose  of  solving  systems  of  linear  algebraic 
equations/  A  method  known  as  the  "Two-step  fraction-free 
Gaussian  elimination"  [14]  has  been  Implemented.  This 
method  Is  an  Improvement  over  a  corresponding  one-step 
method  and  features  a  procedure  that  keeps  the  size  of  the 
Intermediate  expressions  In  the  course  of  the  reduction  down 
by  dividing  them  by  a  common  factor  which  the  procedure  can 
predict.  The  advantage  of  this  method  over  a  more 
efficient  and  elaborate  scheme  [18]  Is  Its  simplicity. 
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6.  Convergence  of  Integrals 

There  are  many  methods  to  determine  the  convergence  or 
divergence  of  a  given  Improper  Integral.  Some  are  limit 
tests#  others  comparison  tests  etc.  Some  test  for  absolute 
convergence#  others  conditional  or  uniform  convergence.  It 
Is  a  possible  area  Tor  future  work.  There  Is  no  real  need 
to  have  such  elaborate  schemes  In  WANDERER#  for  after  all  It 
Is  the  value  of  the  given  Integral  that  Is  desired. 

Consider 


If  f(X)  Is  a  rational  function  of  X#  WANDERER  combines 
convergence  tests  with  evaluation  algorithms  as  explained  In 
chapters  4  and  5.  If  f  Is  not  rational#  then  It  Is  not 
allowed  to  become  Infinite  except  at  a  and  b.  This  Is  a 
convention  on  Inputs  used  by  WANDERER.  Thus#  J  Is  a  proper 
Integral  If  a#  b#  f(a)  and  f(b)  are  finite.  If  J  Is 
Improper#  WANDERER  has  a  test  for  absolute  divergence. 

Test  for  absolute  divergence: 

(1)  If  a  and  b  are  finite  and  b  >  a#  then  J  Is 
absolutely  divergent  If 
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LIMIT((b-X)f(X),X,b, MINUS)  j»  0, 
or 

LIMIT((X-a)f(X)#X#a#PLUS)  t  0. 

(2)  If  a  Is  Infinite#  then  J  Is  absolutely  divergent  If 

UMIT(X*f(X)/X/a>  f  0. 

(3)  If  b  Is  Infinite#  then  J  Is  absolutely  divergent  If 

LIMIT(X*f (X)/X/b)  |»  0. 

If  a  function  G(X)  exists  such  that  G'(X)  *  f(X)/  then 
the  method  of  antiderivative  can  be  used.  That  Is  to  com¬ 
pute  J  by  evaluating 

LIMIT(G(X)/X/b/MIMUS)  -  LIMIT(G(X)/X/a/PLUS). 

If  this  value  Is  finite#  It  Is  the  value  of  the  given 
Improper  Integral.  If  It  does  not  exist#  then  J  diverges. 
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CHAPTER  VI  I 

AN  APPLICATION  OF  MACSYMA 
AND  WANDERER 


0.  Introduction 

The  usefulness  of  a  general  purpose  algebraic  manipula¬ 
tion  system  such  as  MACSYMA  In  facilitating  the  solution  of 
mathematical  problems  has  been  demonstrated.  Using  his 
"Symbolic  Mathematical  Laboratory"  [15],  Martin  demonstrated 
solutions  to  three  demanding  problems  tn  applied 
mathematics.  These  examples  emphasize  the  fact  that  routine 
algebraic  computation  can  be  done  by  computer  programs  not 
only  without  error  but  much  faster  than  by  hand.  The  value 
of  such  a  system  Is  especially  appreciated  when  the  expres¬ 
sions  Involved  are  large  and  complicated.  By  employing 
such  a  computer  facility  the  human  problem  solver  may  be 
freed  from  the  tedious  and  uninspiring  manipulations  to 
think  more  about  the  profound  aspects  of  hls  problem.  But 
this  Is  not  all  such  a  system  can  do.  Moses‘s  SIN  [20],  a 
program  for  Indefinite  Integration,  provides  a  good  example 
of  successful  mechanization  of  a  mathematical  process  which 
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Is  far  from  routine  or  straight-forward.  The  Introduction 
of  SIN  broadened  the  scope  of  algebraic  manipulation  systems 
slgnl f Icantly* 

The  purpose  of  this  chapter  Is  to  show  how  MACSYMA  and 
WANDERER  can  be  used  to  help  solve  complicated  problems  that 
are  of  practical  Importance.  One  such  problem  Is  the  asymp¬ 
totic  evaluation  of  certain  contour  Integrals  arising  In 
mathematical  physics.  Usually  one  starts  with  one  or  a  set 
of  differential  equations  describing  a  physical  problem. 
Solving  these  equations  by  one  method  or  another,  most  often 
by  Integral  transforms,  one  will  arrive  at  a  solution  In  the 
form  of  a  definite  Integral  which  Is  often  difficult.  If  not 
Impossible,  to  evaluate  exactly.  Frequently,  one  Is  not  so 
Interested  In  the  exact  solution  but  the  behavior  of  the 
system  when  one  parameter  becomes  very  large  or  small.  This 
Is  where  asymptotic  analysis  Is  needed* 

MACSYMA  Is  used  to  obtain  the  asymptotic  solution  of  an 
Infinite  integral.  Many  facilities  provided  In  such  an 
algebraic  manipulation  system  can  be  Illustrated  through 
this  application.  It  Is  also  possible  to  show  how  the 
definite  Integration  capability  provided  by  WANDERER  Is 
needed  for  the  successful  computation  of  the  results. 
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1.  An  Asymptotic  Analysis  Problem  and  an  Outline  of  the 
Method  of  Steepest  Descent 


Consider  the  Infinite  Integral 

C  M( t) 

J(C)  -  /  %E  dt. 


rco 

a 

J-ao 


C  >  0 


where 


H(t)  «  -  —  -  %\  t 
4 


While  Integrating  J  exactly  may  be  Impossible,  Its 
asymptotic  behavior  as  C  becomes  very  large  can  be 
Investigated.  To  obtain  the  asymptotic  expansion  of  J,  the 
method  of  steepest  descent  [5,7]  will  be  employed. 

Basically,  the  method  of  steepest  descent  consists  In 
deforming  the  contour  of  Integration  In  such  a  way  that  the 
major  contribution  to  the  Integral  arises  from  a  small  por¬ 
tion  of  the  new  path  of  Integration.  The  contribution  will 
become  more  and  more  dominant  as  the  parameter  of  Interest 
grows.  This  parameter  here  Is  C. 

The  first  step  In  this  method  Is  to  find  the  new  path 
of  Integration.  On  a  given  contour,  larger  contributions 
comes  from  portions  where  the  Integrand  Is  larger  In 
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absolute  value  and  less  oscillatory.  Hence  the  requirements 
for  a  desirable  path  are:  (a)  the  absolute  value  of  the 
Integrand  becomes  maximum  at  a  point/  tl  say,  on  the  path, 
(b)  The  argument  of  the  Integrand  Is  constant  on  the  path 
near  tl.  The  first  requirement  Is  obvious.  The  second  Is 
essential,  for  If  the,  phase  angle  changes  even  slightly  near 
tl,  this  change  will  be  magnified  by  the  very  large  factor  C 
resulting  In  rapid  oscillations  of  the  Integrand  and 
therefore  negattng  any  possible  contribution  from  the  point 
tl.  Let  UCtl,t2)«  Rl(H(tini  t2))  and  V(tl,t2)  »  lm(H(tl+*l 
t2)).  Let  p  be  a  point  where 

dU  dU  dV  dV 

—  —  a  -  -  --  a  -  a  -  -  —  -  a  0* 

d  tl  d  t2  d  tl  d  t2 

Then  p  Is  certainly  a  candidate  for  tl  required  In  (a)  and 
(b).  Such  a  point  Is  called  a  saddle  point.  The  choice  of 
the  name  saddle  point  will  be  made  clearer  later.  The 
Cauchy-R teman  conditions  Imply  that  criterion  for  such  a 
point  Is  H*(t)=0.  There  may  be  more  than  one  such  point  In 
the  complex  t-plane.  There  are  an  Infinite  number  of  diffe¬ 
rent  curves  which  pass  through  a  saddle  point  and  satisfy 
(a)  and  (b).  Among  them  the  path  along  which  the  Integrand 
decrease  In  size  most  rapidly  Is  the  best.  Tor  the  Integral 
J,  this  means  a  curve  on  which  U  decreases  most  rapidly. 
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ri 

) 


Let  H'(b)  “0.  It  can  be  shown  by  use  of  properties  of 
analytic  functions  that  IJ  varies  most  rapidly  on  curves 
lm(H( t) ) °const.  (usually  called  level  curves).  If  H1 
two  level  curves  will  pass  through  the  point  b#  Intersecting 
at  right  angle/  as  shown  In  fig.  1. 


On  one  of  these  two  curves#  say  curve  A#  Rl(H(t))  Is  minimum 
at  b  and  Increases  as  t  moves  along  A  away  from  b.  On  the 
other  hand#  Rl CH(t>>  Is  maximum  at  t=b  on  the  curve  D  and 
decreases  as  t  moves  away  from  b  on  D.  Curve  A  Is  called 
the  steepest  ascent  path  and  D  the  steepest  descent  path. 

The  point  b  Is  usually  referred  to  as  a  saddle  point.  If 
HM(b)j*0#  b  Is  a  saddle  point  of  order  1#  If  HM(b)B0  and 
HM,(b)j*0  of  order  2#  etc.  If  the  original  contour  can  be 
deformed  onto  one  or  a  combination  of  such  steepest  descent 
paths#  then  the  asymptotic  expansion  of  the  given  Integral 
can  be  obtained  by  a  rather  routine  procedure  which 
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Involves:  change  of  variable  of  Integration#  Inversion  of 
truncated  power  series  at  the  saddle  points  and  term  by  term 
Integration. 

Let  us  outline  the  steps  of  the  solution  procedure  as 
fol lows. 

1)  Locate  And  determine  the  order  of  saddle  points 
of  H(T) . 

2)  Compute  U(tl#t2)  and  V(tl,t2)  such  that 
H(tl  +  %\  t2)  =  U(tl,t2)  ♦  %l  V( tl# t2 ) , 

3)  Obtain  V(tl,t2)  3  const,  curves  .which  pass 
through  the  relevant  saddle  points. 

4)  Examine  the  V(tl/t2)  *  const,  curves  to 
determine  whether  deformation  of  contour  can  be  made 
to  curves  through  the  saddle  points. 

5)  Change  the  variable  of  Integration. 

6)  Express  t  as  a  truncated  series  In  the  new 
variable  about  each  relevant  saddle  point. 

7)  Determine  the  coefficients  In  the  above  series. 

8)  Apply  Watson's  lemma  to  obtain  the  first  few 
terms  of  the  asymptotic  expansion  by  Integrating 
term  by  term. 
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2.  Solution  Steps 

Presented  here  are  the  solution  steps  of  the  asymptotic 
expansion  of  J(C)  In  the  exact  sequence  as  they  have  been 
carried  out  using  MACSYMA .  The  lines  labelled  (Cl)  are 
Input  commands.  A  command  line  ends  with  either  a  Q  sign  or 
a  $  sign.  The  Q  sign  causes  results  obtained  by  executing 
the  command  line  to  be  displayed  In  a  subsequent  line 
labelled  (Dl).  A  (Cl)  line  together  with  a  corresponding  CVi) 
line  will  be  referred  to  as  step  I.  The  $  as  an  end  of 
command  line  character  suppresses  display  of  results  for 
that  line.  Explanatory  texts  will  be  Inserted  between 
lines.  To  avoid  becoming  a  user's  manual  for  MACSYMA, 
explanation  for  the  commands  used  In  the  solution  will  be 
made  quite  brief.  For  a  more  detailed  look  at  MACSYMA  the 
reader  Is  referred  to  [16], 

(Cl)  P: 8$ 

P  Is  a  parameter  which  Is  set  depending  on  the  number 
of  terms  desired  In  the  asymptotic  expansion.  By  setting  P 
to  8,  we  shall  obtain  the  first  4  terms.  The  reason  will 
become  evident  later. 

(C2)  H(T) :a-T**4/4-%l  *T@ 
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(02) 


H(T> : 1 


defining  the  functton  H(T) 


T 


(C3 )  DIFF(H(T),T)=0@ 
(D3) 


3 

-  T  -?l»0 


creating  an  equation  H'(T)»0 


(C4 )  S0LVE(5,T)Q 
SOLUTION 

SQRT( 3 )  -  51 

(E4)  T  -  . 

2 

(E5)  T  «  51 

-  SQRT(3)  -  51 

(E6)  T  =  . 

2 

(D6)  ( E4, E5, T6) 

The  5  sign  used  In  (C4)  stands  for  the  last  (03),  In 
general  a  %  sign  represents  the  last  expression  labelled 
(01).  SOLVE  In  (C4)  Is  an  Invocation  of  the  MACSYMA  SOLVE 
program  (see  [10]  for  Its  capabilities  and  limitations). 
The  roots  of  H'(T)a0  give  three  first  order  saddle  points. 
Passing  through  each  of  these  points  there  will  be  one 
steepest  ascent  and  one  steepest  descent  level  curve 
(Im  (H(t))»  const,  curves).  The  saddle  points  in  (E4)  and 
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(E6)  are  the  points  B  and  A  respectively  shown  In  fig.  2. 

The  value  of  H  at  these  three  saddle  points  will  now  be 
computed. 

( C 7 )  HB:RATSIMP(H( PARK E4,2)))@ 

-  3  %\  SQRTC3)  -  3 

(D7)  . 

8 

What  has  been  done  In  step  7  Is  the  computation  and 
simplification  of  H(B).  The  command  RATS  IMP  causes  rational 
simplification  which  Is  essentially  putting  expressions  to 
be  simplified  Into  the  form  of  one  numerator. and  one 
denominator  and  perform  all  possible  GCD  cancellations.  The 
command  PART  allows  a  user  to  obtain  subexpressions  of  an 
expression.  PART(E4/2)  returns  the  second  part  of  the 
equation  E4  which  Is  Its  right-hand  side.  Mote  that 
commands  can  be  nested. 


(C8)  HC:RATSIMP(H(PART(E5,2)))@ 


(08) 

(C9)  HA: RATS IMP(H( PART( E6, 2) ) )@ 


3 

4 


(D9) 


3  %\  SQRT( 3)  -  3 
8 


The  point  T-%\  turns  out  to  be  Irrelevant  because  the 
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new  path  which  will  be  determined  does  not  pass  through  It. 
The  next  goal  Is  to  obtain  curves  passing  through  the  saddle 
points  (SQRT(3)/2,  -1/2)  and  (-SQRT(3)/2,  -1/2)  along  which 
the  Imaginary  part  of  H(T)  Is  constant  and  the  real  part  of 
H(T)  varies  most  rapidly,  l.e.  the  steepest  paths  through 
the  saddle  points. 

(CIO)  EXP: EXPAND (H(Tl+§ I *T2)  )@ 

4  2  2 

T2  3  3  T1  T2 

(DIO) - ♦  %\  T1  T2  + . - 

4  2 

3 

-§l  T1  T2  ♦  T2 

(Cll)  V(T1,T2):«"COEFF(D1O,§I)0 

3  3 

(Dll)  V(T1,T2): 3  T1  T2  -  T1  T2  -  T1 

By  steps  10  and  11,  the  Imaginary  part  of  H(T1+  %\  T2) 
Is  found  and  given  a  function  name  V(T1,T2).  The  command 
C0EFF(exp, var,n)  computes  the  coefficient  of  var**n  In  exp. 
In  this  case,  the  coefficient  of  %\  In  (DIO)  Is  exactly  the 
Imaginary  part  of  (DIO). 


T1 

. %\  T1 

4 


(C12)  U(T1,T2):=,,C0EFF(EXP,  §1,0)0 
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4  2  2  4 

T2  3  T1  T2  T1 

(D12)  U(T1,T2):  = - ♦ . +  T2 - 

4  2  4 

U(T1/T2)  Is  defined  to  be  the  real  part  of 

H(T1  ♦  *1  T2>. 


(C13)V(Tl,T2)=V(SQRT(3)/2,  -1/2)0 

33  3  SQRT( 3) 

(D13)  T1  T2  -  T1  T2  -  T1  •  . 

8 

Obtained  In  (D13)  Is  the  equation  of  a  curve,  V  ■ 
constant, which  passes  through  the  saddle  point  B, 
(SQRT(3)/2,  -1/2).  This  curve  will  be  referred  to  as  CR. 

(C14)V(Tl,T2)=V(-SQRT(3)/2/  -1/2)0 

33  3  SQRT( 3) 

(D14)  T1  T2  -  T1  T2  -  T1 . 

8 

(014)  Is  the  equation  of  a  curve,  V=constant,  passing 
through  the  saddle  point  A,  (-SQRT(3)/2,  -1/2).  Let  this 
curve  be  CL.  The  curves  CL  and  CR  have  to  be  examined 
carefully  by  the  human  problem  solver  to  determlng  the  new 
path  of  Integration.  The  manner  In  which  they  extend  to 
Infinity  Is  often  Important  In  deforming  the  contour. 
Asymptotes  to  these  curves  can  be  found  easily.  For  an 
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algebraic  curve#  f(X#Y)aO#  the  way  to  find  asymptotes  Is  to 
substitute  m*X+c  for  Y  In  f(X,Y)»0  then  determine  values  of 
m  and  c  such  that  the  equation  has  two  Infinite  roots  (l.e.# 
equating  to  zero  the  coefficient  of  the  highest  and  second 
highest  powers  of  X).  Asymptotes  parallel  to  X«0  are  missed 
by  this  method.  f(X#Y)a0  has  such  an  asymptote  If  It  Is 
possible  to  choose  h  In  such  a  way  that  the  equation 
f(h#Y)a0  has  two  Infinite  roots.  Asymptotes  will  be  found 
for  CL  and  CR  in  the  next  two  steps. 

(C15)  SUBST I TUTE( (T2aA*Tl+B ) # PART (D 14# 1) ) 0 

3  3 

(015)  -  T1  (A  T1  ♦  B)  ♦  T1  (A  T1  ♦  B)  -  T1 

(CIO  RATSIMP(%)Q 

3  4  2  3  2  2  3 

(D16)(A  -  A)  T1  ♦  (3  A  -  l)B  T1  ♦  3  A  B  T1  ♦  (B  -  1)  T1 

The  asymptotes  are  clearly  Tla0#  T2=0#  T1=T2  and  Tl*  - 
T2.  With  the  aid  of  MACSYMA  to  generate  points  we  plotted 
the  curves  as  shown  In  fig.  2.  with  steepest  descent  paths 
labelled  Ca  and  Cb.  This  figure  Is  not  a  computer  generated 
plot. 


*L-tZ 


-n 
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It  can  be  seen  In  fig.  2  that  by  use  of  Cauchy's 
Integral  theorem/  the  original  contour  can  be  deformed  to 
the  contour  Ca+Cb.  For  the  purpose  of  change  of  path  of 
Integration  a  new  variable  R  Is  Introduced  which  will  be  the 
parameter  of  our  steepest  descent  paths. 


(C17)  -R**2=H(T)-H(T0)Q 

2 

(017)  -  R  «  V  TO 


4 

T 

4 


T  ♦ 


1 

4 


TO  stands  for  either  of  the  two  saddle  points  A  or  B. 
Indeed  If  T0«A/  R  Is  real  tf  and  only  If  T  Is  a  point  on  Ca/ 
for  lm(H(T)-H( A) )a0  and  0  S  (H(T)-H(A))  only  for  points  on 
0*.  The  same  can  be  said  about  T0»B.  T  has  to  be  expressed 
In  terms  of  R  In  order  to  perform  the  change  of  variable. 
Solving  for  T  as  a  function  of  R  exactly  In  (D17)  Is  not 
necessary.  What  Is  needed  Is  the  first  few  terms  of  a  power 
series  expansion  of  T  In  terms  of  R  about  the  point  T**T0 
(the  saddle  point). 


(C18)T=TO+DOSUM(J/1,P/B[J]*R**J)@ 


8  7  6  5  4 

(D 18 )T  3  TO  ♦  B  R  +BR  ♦  B  R  ♦  B  R+B  R 
8  7  6  5  4 

3  2 

♦  B  R  +B^R  ♦  B,  R 
3  2  1 
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T  Is  set  to  a  truncated  power  series  In  R.  The  unknown 
coefficients  BI  will  be  solved  for  by  use  of  (D17). 

(C19)  E:SUBSTITUTE(D18,  R**2+PART(D17,2))@ 

8  7  6  5  4 

(D19)  -  (TO  ♦  B  R  ♦  B  R  ♦  B  R  ♦  B  R  +  B  R 

8  7  6  5  4 

3  2  4  8  7 

♦  B  R  +  B  R  ♦  B  R)  /4  -  %\  (TO  ♦  B  R  ♦  B  R 

3  2  1  8  7 

6  5  4  3  2 

+  BR+B  R+BR+BR+BR+BR) 

6  5  4  3  2  1 

2  1 

+  V  TO  ♦  R  ♦  - 

4 

The  expression  (D19)  has  been  given  a  name  E.  Before  E 
Is  expanded  In  order  to  collect  terms,  a  few  simplification 
rules  are  defined  so  as  to  discard  powers  of  R  higher  than  8 
In  the  expansion  process.  This  greatly  reduces  the 
Intermediate  expression  bulge  which  would  otherwise  occur. 


(C20)  DECLARE  <N,NPRED)$ 

(C21)  MPRED(X) : a I F  X  >  P  THEN  TRUE  ELSE  FALSES 

(C22)  TELLSIMP  (R**N,0)$ 

0 

IS  THE  REPLACEMENT 


A  simplification  rule  has  been  set  up  so  that  any  R**n 
Is  replaced  by  0  If  N  >  p  which  Is  8. 
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(C23)  TELLSIMP(T0**3/-%I)$ 

An  additional  simplification  rule  Is  defined  so  that 
T0**3  shall  be  replaced  by  -%l  automatically. 

(C24 )  RATVARS( R) $ 

R  shall  be  the  main  variable  In  subsequent  rational 
slmpl If Icatlons. 

(C25)  E:  PARK  RAT(  E) / 1)0 
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(025) 

8  2  2 
R  ((  -  12  B  B  -  12  B  B  -  12  B  B  -  6  B  )  TO 
1  7  2  6  3  5  4 


♦  (  -  12  B  B  -  24  B  B  B  ♦  (  -  24  B  B  -12 
16  12  5  13 

2  3  2  2 

■  12  B  B  )  TO  •  4  B  B  -  12  B  BB-6B 
2  3  1  5  1  2  4  1 


2  4 

-  12  B  B  B  -  B  ) 

12  3  2 

7  2 

♦  R  ((-12  B  B  -  12  B  B  -12B  B)TO 
1  6  2  5  3  4 

2  2  2 
♦  (  -  12  B  B  -  24  B  B  B  -  12  B  B  -  12  B 


1  5 


12  4 


1  3 


3  2  3 

-  4  B  B  -  12  B  BB-4BB) 
14  12  3  12 

6  2  2 
♦  R  ((-12  B  B  -12B  B  -6B  )T0 

1  5  2  4  3 

2  3 

+  (  -  12  B  B  -24B  B  B  -  4B  )T0 
14  12  3  2 

3  2  2 

-48  B-6B  B) 

13  12 


5  2 

+  R  ((  -  12  B  B  -  12  B  B  )  TO 
14  2  3 


2 

B  )  B 

2  4 

2 

B 

3 


B  )  TO 
3 
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2  2  3 

♦  (  -  12  B  B  -  12  B  B  )T0  -4B  B) 

13  12  12 


4  2  2  2  4 

+  R  ((  -  12  B  B  -  6  B  )  TO  -12B  B  TO  -  B  ) 
13  2  12  1 

3  2  3  2  2  2 

♦  R  (  -  12  B  B  TO  -  4  B  TO)  +  R  (4  -  6  B  TO  ) 
1  '2  1  1 


E  Is  set  to  the  huge  expression  above  which  Is  the 
numerator  of  a  truncated  expansion  of  E;  the  denominator  Is 
a  constant. 


(C26)F0R  J: 1  STEP  1  UNTIL  d  >  (P-1) 

DO  BLOCK  (EQ[J] :COEFF(E,R, J+1)=0,DISPLAY(EQ[J] ))$ 

This  do  loop  Is  used  to  generated  the  coefficients  of 
the  various  powers  of  R  In  E  and  store  them  In  the  form  of 
equations  In  an  array  RL.  Each  EQ[J]  will  he  displayed 
after  It  Is  set. 


2 

4  -  6  B 

1 

2 

-  12  B  B  TO 
1  2 


2 

TO  =  0 


3 

-  4  B  TO  »  0 
1 


2  2  2  4 

-  6  3  )  TO  -  12  B  B  TO  -  B 

2  12  1 


0 


(  -  12  B  B 
1  3 
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(  - 

♦ 

(  - 

♦  (  -  12 

(  - 

♦  (  -  12 

B  -  12 
4 


12  B  B  -  12  B  B  )  TO 
14  2  3 


♦  (  -  12  B  B  -12B  B  )  TO  -  4  B  B  »  0 


1  3 


1  2 


1  2 


2  2 

12  B  B  -  12  B  B  -  6  B  )  TO 
1  5  2  4  3 


C-12B  B  •  24  B  B  B  -4B  )TO 
14  12  3  2 

3  2  2 

•  4  B  B-6B  B  ■  0 
13  12 


12  B  B  -  12  B  B  -  12  B  B  )  TO 

1  6  2  5  3  4 

2  2  2 
B  B-24BBB-12BB  -  12  B  B) 
1  5  1  2  4  1  3  2  3 

3  2  3 

TO-4B  B  -  12  B  BB-4BB  »  0 
14  12  3  12 


2  2 

12  B  B  -  12  B  B  -  12  B  B  -6B  )  TO 

1  7  2  6  3  5  4 

2  2 
B  B  -  24  B  B  B  +  C-24B  B  -12B  ) 

16  125  13  2 

2  3  2 

B  B  )  TO  -  4  B  B  -  12  B  BB 
2  3  15  12  4 
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2  2  2  4 

-6B  B  -  12  B  B  B-B  »  0 
13  12  3  2 

Being  conscious  of  the  lack  of  storage  space  left  for 
subsequent  man Ipulat Ions,  we  get  rid  of  the  computation 
history  to  create  some  space. 

(C27)  K I LL (H I  STORY) @ 

(Cl)  S0LVE(EQ[1J,B[1])@ 

SQRT(6) 

(El)  B  . 

1  3  TO 

SQRT( 6)  • 

(E2)  B  . . 

1  3  TO 

(02)  (El, E2) 

The  choice  of  the  value  for  B1  here  affects  only  the 
sense  with  which  the  steepest  decent  paths  are  traversed. 
Thus  If  we  choose  (El),  R  would  vary  from  IMF  to  MINF  on  Ca 
and  from  MINF  to  INF  on  Cb. 

( C 3 )  B [1] : PART (El, 2)$ 

A  do  loop  Is  used  In  (C4)  to  solve  for  the  remaining 
B's  (B2  through  B7). 

(C4)  FOR  J: 2  STEP  1  UNTIL  J  >  (P-1) 

DO  BLOCK  (SOL:  EV (SOLVE ( EQ  [J]  ,B[J]  ),EVAL), 

B[JJ :PART($0L,2),DISPLAY(S0L))$ 

The  EV  command,  with  EVAL  as  the  third  argument,  causes 
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the  answer  returned  by  SOLVE  to  be  evaluated  once  more. 
Since  every  BJ  Is  solved  In  terms  of  earlier  B's  which 
already  have  some  value/  one  more  level  of  evaluation  will 
eliminate  the  dependencies  on  previous  B's, 

2  %\ 

B  B  *  “““* 

2  9 

10  %\  -  3  TO 

B  . . 

3  2 

27  SQRT( 6)  TO 


9  %\  TO  +  23 

B  - - 

4  243  TO 

2 

-  3  %\  TO  ♦  28  TO  -  108  %\ 

B  . 

5  972  SQRT(6) 


B 

7 


2 

225  %\  TO  -  342  TO  ♦  1031  %\ 

B  . 

6  2 
26244  TO 

2  2 

5  %\  TO  SQRT( 6)  %\  TO  47  SQRT(6)  TO 

-  ------------  +  --------------  +  ------------- 

1458  SQRT(6)  729  6561 

107  TO  6839  %\  170  SQRT(6)  %\ 

+  -----------  - - -  -------------- 

2187  SQRTC6)  52488  SQRT(6)  19683 


1927  %\ 
26244  SQRT( 6) 
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The  unknown  coefficients  Bl  have  been  determined.  The 

next  goal  Is  to  compute 

®  2  dT 

EXPC-C  R  )  —  dR 


/• 


dR 


dT 

for  T0*A  and  T0=B,  Odd  terms  of  R  In  --  do  not  contribute, 

dR 

Therefore  only  4  terms  need  be  Integrated.  These  terms  are 
even  In  R,  thus  the  range  of  Integration  can  be  changed  to 
(0  IMF).  A  function  FN(X)  will  be  defined  for  carrying  out 
this  term-by-term  Integration  at  a  variable  point  TOaX, 


(C5)  FM(X):=BL0CK(AN5:C/  FOR  J:1  STEP  2  UNTIL  J  >  (P-1) 

DO  ANS: AMS +J* RATS IMP( EV(B | J | ,TOaX) )* 

DEFINTUE**(-OR**2)*R**(  J-1),R,0,  INF)  ,  RETURN  CANS))  $ 

(CG)  FNC-SQRTC 3)/2-%l /2)@ 

IS  THE  EXPRESSION 
-  C 

POSITIVE/  NEGATIVE/  OR  ZERO 
NEGATIVES 

For  the  purpose  of  Integration/  WANDERER  askes  the  user 
about  the  sign  of  C.  The  contribution  from  saddle  point  A 
Is  given  In  the  next  result. 


(DG) 


SQRTC6)  SQRTC  %PI ) 

(3  SQRTC  3)  +  3  %\ )  SQRT(C) 


I 
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3  SQRTUPI)  (3  SQRT( 3)  +  23  $1) 


3/2 

4  SQRT( 6)  (27  %\  SQRT(3)  ♦  27)  C 

5  SQRT( %PI )  (25  SQRT( 3)  ♦  247  ^ I ) 


5/2 

5184  SQRT(6)  C 

-  (35  SQRT(J?;PI)  (SQRT(6>  (672  SQRT(3)  ♦  1816  %\) 

3582  SQRT(  3)  13077  %\  7/2 

. . ))/( 839808  C  ) 

SQRT( 6)  SQRT( 6) 


(C7)  FN(SQRT(3)/2-*l/2)@ 

IS  THE  EXPRESSION 
-  C 

POSITIVE,  NEGATIVE,  OR  ZERO 
NEGATIVES 


SQRT(  6)  SQRT(fcPI) 

(D7)  . 

(3  SQRT(3)  -  3  %l)  SQRT(C) 

3  SQRTUPI)  (3  SQRT( 3)  -  23  *1) 

+  -------------- — ------------------ 

3/2 

4  SQRT( 6)  (27  %\  SQRT(3)  -  27)  C 

5  SQRTUPI)  (25  SQRT( 3)  -  247  %\ ) 

+  -------------- — - —  - - 

5/2 

5184  SQRT(6)  C 


♦  (35  SQRT(|PI)  (SQRT( 6)  (672  SQRT(3)  -  1816  %\) 

3582  SQRT( 3)  13077  %\  7/2 

♦ . )  )/( 830808  C  ) 

SQRT(G)  SQRT( 6 ) 
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The  above  Is  the  contribution  from  saddle  point  B.  Now 
the  final  result  Is  obtained  by  computing  the  following 

(-3  yr  %i -3)/8  (3  yr  ?,i-3)/8 

2  UE  D7  -  2  *E  D6) 


where  D6  and  07  are  as  given  above.  This  expression  Is,  by 
Inspection,  equivalent  to 

(-3  yr  % i -3)/8  (3  yr  *i-3)/8  _ 

2  UE  07  ♦  %E  07) 

where  D7  Is  the  complex  conjugate  of  07.  Therefore, 


(C8)  4*RL( %E**HA*07)@ 

3  SQRT( 3)  %\  -  3 

8  SQRTC6)  SQRTUPI) 

( . 

(3  SQRT( 3)  -  3  *1)  SQRT(C) 

3  SQRTUP I )  (3  SQRT( 3)  -  23  *1) 


(08)  4  RL(%E 


3/2 

4  SQRT( 6)  (27  %\  SQRT(3)  -  27)  C. 

5  SQRTUPI)  (25  SQRT(  3)  -  247  %\) 


5/2 

5184  SQRT(6)  C 


+  (35  SQRTUPI)  (SQRT( 6)  (672  SQRT(3)  -  1816  %\) 


3582  SQRT( 3) 


SQRT( 6) 


13077  %\  7/2 

. ))/( 839808  C  )) 

SQRT(6) 
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This  Is  the  first  four  terms  of  the  desired  asymptotic 
expansion,  RL  Is  not  a  command  of  MACSYMA.  It  Is  used  here 
to  denote  the  real  part  of  an  expression. 
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CHAPTER  VIII 

CONCLUSIONS  0  SUGGESTIONS  FOR  FURTHER  WORK 

The  evaluation  of  definite  Integrals  Is  a  classical 
problem  In  mathematics.  Great  Ingenuity  Is  frequently 
requlred/  with  many  Integrals  demanding  special  devices. 

The  lack  of  a  sufficiently  general  theory  makes  evaluation 
of  definite  Integrals  very  difficult.  It  Is  doubtful  that  a 
theory#  comparable  In  generality  to  the  Rlsch  Integration 
algorithm  [26]  for  Indefinite  Integrals#  can  be  developed  In 
the  near  future.  We  have  shown  that  the  convergence  of  a 
class  of  Integrals  of  elementary  functions  Is  recursively 
undecldable.  The  proof  Is  In  appendix  E.  The  WANDERER  pro¬ 
gram  presented  here  Is  a  prototype  heuristic  computer  pro¬ 
gram  for  the  symbolic  evaluation  of  definite  Integrals.  It 
contains  both  general  methods  such  as  contour  Integration# 
residue  theory  and  differentiation  with  respect  to  a 
parameter#  In  addition  to  quite  a  few  special  methods  for 
specific  types  of  Integrals.  Clues  as  to  which  method  to 
use  for  a  given  problem  are  obtained  from  the  Integration 
range  and  the  form  of  the  Integrand.  Although  many  types  of 
definite  Integrals  can  be  evaluated  by  WANDERER#  It/  as 


•  ■  j* '  h*  iiv-'i  SV?  v  ^vt  *  lyf  ^Jfr^'iV.*  -  •'^V* 


m  **  •>  c  r*‘iP!»  iMUpttTr  “ 
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almost  any  other  heuristic  computer  program,  has  Its  limita¬ 
tions.  Foremost  among  these  Is  the  fact  that  It  Is  a  pro¬ 
gram  for  the  evaluation  of  real  definite  Integrals  of 
elementary  functions. 

It  Is  hoped  that  the  work  reported  here  may  provide  a 
starting  point  for  new  approaches  to  the  evaluation  of 
definite  Integrals  from  the  viewpoint  of  symbol  manipula¬ 
tion.  The  advantage  of  this  approach  Is  twofold:  (1)  A  com¬ 
puter  can  use  Integration  methods  that  are  too  lengthy  or 
complicated  to  be  carried  out  by  hand;  (2)  Such  a  computer 
program  contains  a  collection  of  powerful  methods  that  can 
Interact  with  one  another  and  can  produce  answers  to  Inte¬ 
grals  not  present  In  any  finite  table.  For  this  reason#  a 
good  definite  Integration  program  together  with  other 
facilities  provided  In  an  algebraic  manipulation  system  can 
be  very  useful  to  applied  mathematicians. 

DELIMITER  Is  a  rather  sophisticated  program  for  comput¬ 
ing  limits  of  elementary  functions.  Such  a  program  has  been 
shown  to  be  a  very  useful  tool  in  an  algebraic  manipulation 
system.  The  method  of  comparing  orders  of  Infinity  used  In 
DELIMITER  Is  an  Important  concept  which  Is  useful  In  places 
other  than  the  computation  of  limits. 
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A  basic  assumption  of  DELIMITER  Is  that  the  operator 
LIMIT  Is  distributive  over  the  operators  ♦,  -,  *,  and  **, 
l.e .,  rules  such  as  the  limit  of  a  sum  Is  the  sum  of  the 
limits  hold.  These  rules  are  valid  only  If  the  subproblems 
thus  generated  produce  answers  which  do  not  lead  to  an 
Indeterminate  form.  Of  course  there  Is  L'Hospltal's  rule 
which  can  be  applied  to  some  Indeterminate  forms,  but 
generally  only  those  Involving  only  the  operator  *  or  **. 
Although  many  aspects  of  the  Indeterminate  form  (INF-INF) 
have  been  considered.  In  chapter  2,  the  algorithm  for  Its 
determination  Is  not  complete.  A  powerful  series  expansion 
program  would  be  helpful  In  some  cases.  Yet,  such  a  pro¬ 
gram  can  not  solve  all  problems.  Consider  for  Instance,  the 
limit  problem 


2  2 

L I M I T( S I N  X  ♦  COS  X,X, INF) . 

To  obtain  such  limits,  an  algebraic  manipulation  system 
must  be  able  to  detect  all  constant  Identities.  This  Is  not 
possible  for  the  set  of  all  expressions  (2  2] .  It  may  be 
possible  for  a  proper  subset  of  all  expressions.  For 
example,  many  trigonometrical  Identities  disappear  If  all 
trigonometric  functions  ore  radically  transformed  Into  sums 
of  complex  exponentials.  This  Is  a  basic  problem  of  great 
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practical  and  theoretical  Interest  [6],  Such  simplification 
capabilities  would  be  of  use  to  many  other  programs  as  well. 

DELIMITER  can  be  augmented  by  Increasing  the  types  of 
functions  It  can  handle  which  may  Include  functions  defined 
by  Integrals.  Another  possible  area  of  research  Is  the 
automatic  determination  of  superior  and  Inferior  limits. 

The  fact  that  WANDERER  can  evaluate  many  non-trlvlal 
definite  Integrals  does  not  mean  that  It  can  compete  with  an 
expert  human  Integrator  yet.  For  one  thing,  a  mathematician 
can  usually  construct  a  function  of  a  complex  variable  and  a 
suitable  contour  for  evaluating  different  Integrals. 

WANDERER  cannot  form  a  contour  based  on  analysts  of  a  given 
real  Integral.  It  simply  selects  from  the  cases  known  to 
It.  WANDERER  would  be  much  more  powerful  If  It  could 
determine,  for  a  given  Integral,  whether  the  method  of  con¬ 
tour  Integration  and  residue  theory  were  applicable,  and.  If 
It  were,  evaluate  the  Integral  by  forming  a  suitable  Inte¬ 
gral  around  a  closed  contour. 

The  evaluation  of  contour  Integrals  by  residue  theory 
usually  requires  the  solution  of  algebraic  or  transcendental 
equations.  WANDERER  uses  the  SOLVE  program  of  MACSYMA  for 
this  purpose.  SOLVE  has  Its  limitations  and  Its  Improvement 


<’<i  <\  ■>' 
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Is  also  a  possible  area  for  further  work. 

Suppose  one  wants  to  evaluate  a  closed  contour  Integral 
by  residue  theory.  An  Interesting  general  question  Is:  what 
knowledge  about  the  Integrand  or  the  functions  used  In  form¬ 
ing  It  Is  necessary.  We  think  the  following  are  essential: 

(1)  evaluation; 

(2)  differentiability  and  derivatives; 

(3)  singularities; 

(4)  asymptotic  behaviour. 

A  natural  extension  of  the  work  here  Is  Integration 
over  arbitrary  user-specified  contours.  This  should  not  be 
difficult  to  do,  except  for  the  lack  of  notation. 

Specif  leal  1  y,  we  must  specify  an  arbitrary  contour  to  a  com¬ 
puter  and  devise  a  general  data  structure  for  use  In 
representing  contours.  For  Instance#  consider  the 
specification  of  the  following  Indented  contour. 


*  * 


v#*  ,-ysy.  WJ*.  ’  y'y.  .; *1?;}  .*/«*»’* **• 
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One  obvious  way  Is  to  specify  a  contour  piecewise. 

Each  piece  of  curve  has  a  parametric  form  and  a  starting  and 
ending  value  for  the  parameter.  There  may  be  other 
approaches. 

Further  work  In  this  area  of  symbol  manipulation  might 
Include: 

(1)  Design  of  computer  algorithms  for  testing  convergence 
and  divergence  of  Integrals, 

(2)  Summation  of  Infinite  series  by  residue  computations, 

(3)  Investigation  of  algorithms  for  definite.  Integration  of 
special  functions. 

(4)  Evaluation  of  multiple  definite  Integrals. 

(5)  Application  of  definite  Integration  programs  In  solution 
of  differential  equations. 
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APPENDIX  A 

A  FLOWCHART  OF  DELIMITER 

A  flowchart  Is  presented  In  the  next  page  which  details 
the  flow  of  control  of  DELIMITER.  Listings  of  the  programs# 
written  In  the  LIS?  programming  language#  may  be  obtained 
from  the  author. 

The  routine  LIM#  appearing  In  the  flowchart#  Is  a 
program  which  applies  the  ’trivial'  limit  rules#  makes  use 
of  subroutines  to  compute  limits  of  sums#  products#  powers# 
and  the  functions  SIN#  COS#  TAN#  LOG#  SINH#  COSH#  TANH.  It 
calls  LIMIT  recursively  and  makes  use  of  L'Hospltal's  rule 
and  other  routines  when  needed. 
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LIMIT(E(X),X,L) 


l'hospital's  rule 

FAI LED?  _ 


E,  X,  U 

.  ~T  , 

E  INDEPENDENT! 

? 

Ireturn  El 


E  RATIONAL 

U LUL 


RATIONAL 

FUNCTION 

ROUTINE 


i — IL *  HOSPITAL'S  RULE 


IRETURN 

COMPARE  ORDERS 
OF  INFINITY  OF 
N  AND  D 

N  AND  D  CONTAIN 
FUNCTIONS  IN  THE 
FORM  e**R(X) 
WHICH  TEND  TO  INF? 


Ur 


RETURN  UM(  NL/DL )l 


N-NUMERATOR 
OF  E 

D-DENOMI NATORl 
■Of  L 


NL  AMD  DL=  /ft/PIVITV?! 

,  tf 

NL  = I NF  OR  MINF  * 
AND  ?  " 

DL=I NF  OR  MINF 

~£~ 

[dlH?F-® 


flP-»  Q  +  ?l 


fE?0 


NL=0  AMD  DL=0fl 


N  AND  D  ARE 
RP-FORMS?  _ 

T 

RP-FORM 

ROUTINE 

ID1 


1HL>0 


IRETURN  INFl 


Ireturn  hinfI- 


RETURN  INFINITYI 


NLH.IMI  T(N,X,  L) 
DL«-LIMI  T(D/X,L) 
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APPENDIX  B 
A  PROOF 


Infinite  Integrals  of 

n 

%\  k  Z  m 
f(Z)  «  %E  Z 

where  n  >  0,  k  i*  0  real,  Rl(m)  >  -1,  and  n  -  Rl(m)  >  1, 
have  been  dlscused  In  chapter  4-3.2.  Some  results  derived 
there  depends  on  the  proof  of  (1)  and  (2)  below. 


The  objective  here  Is  to  supply  the  proof  of 


and 


1)  1 Imlt  f  f(Z)  dZ  «  0, 

o*  J<Lt 

2)  limit  f  f(Z)  dZ  »  0. 

Jen 


where  CR  and  C  are  the  circular  contours  given  below. 
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Let 


le 

Z  «*  r  e  and  m  ■  a  ♦  lb. 


then 


I L 


f(Z)  dz 


(i) 


_  rn/zn  , 

m+l  /  n  lue 

I  r  J  EXP(  I  kr  e  +  i  (m  +  i)  e  )  de 

a+1  rTT/211  n 

=  r  /  EXP(  -  k  SlN(n  e)  r  -  b  e)  de 


M 


For  0  <  e  ^  %P I / C 2  n)  we  have 


EXP(  -  k  S t N(n  e)  r  )  ^  1. 


Thus 

M  s  r 

Since  a  +  1  >  0/  this  completes  the  proof  of  (1). 


,  r'H/zn 
a+1 /  -be 

Jo  *  d* 


Now  from  equation  (l)  we  have 
a+1 


M  = 


'  *  r'K/Z 

n  J0 


n  be 

EX P (  -  k  SIN(e)  r - )  de 

n 


a+1 


2  k  e  r  be 

EX P( - )  de 

fcPI  n 
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$ 


a+1 

%P  I  r 

- (1  -  EXP( 

n 

2  k  n  i*  ♦  %P\  b 


n 

k  r 


The  fact  that  n-a-1  >  0  completes  the  proof  of  (2) 


SPI  b 

- )) 

2  n 
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APPENDIX  C 

PERFORMANCE  OF  WANDERER 
AND  DELIMITER 

As  a  measure  of  the  performance  of  WANDERER  and  the 
limit  programs#  selected  problems  have  been  timed.  The  time 
sharlr.,  system  used  Is  the  ITS  of  the  Artificial 
Intelligence  Laboratory  at  MIT  which  uses  a  PDP-10  computer 
with  a  memory  cycle  time  of  about  2.75  microseconds.  Time 
used  for  parsing  the  Input  string  and  display  of  the 
computed  result  has  been  excluded  In  order  to  obtain  an 
approximation  to  the  time  actually  spent  Inside  WANDERER  or 
DELIMITER.  Garbage  collection  In  the  LISP  [17]  system  In 
which  MACSYMA  Is  written  may  take  place  during  a 
computation.  Although  garbage  collection  Is  a  slow 
procedure#  It  Is  only  fair  to  regard  It  as  part  of  the 
computation  process  being  carried  out.  Therefore#  the  time 
It  required  has  been  Included  In  the  timing  experiments. 

The  results  are  put  In  the  form  of  two  tables.  An  *  Is  used 
In  the  tables  to  Indicate  computation  requiring  LISP  garbage 
col  lection. 
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Limit 


LIMIT 
X  — 0  + 


LIMIT 

X-*+e» 


LIMIT 

X-*+« 


LIMIT 

X-*+oo 


LIMIT 
X-*-  2- 


LIMIT 

X-ff/2 


TABLE  I 

PERFORMANCE  OF  DELIMITER 


Problem  computed 


Resul t 

equivalent  to 


X 

X  1 


X 

(1+1/X) 


X 

X  LOG(X)  +  %E 


3 

4  SQRTCX  +  1) 

LOG( X  +  X  ♦  1)  ♦  %E 


0 


2  2 
X  SQRT(X  +1)  X 

*E  -  %l  IMF 


1 

3  2 

X  -  G  X  ♦  11  X  -  6 
1  -  SIN(X) 

COS(X) 


INF 


0 


Time  t n 
Sec. 

0.82 

1.42 

3.00 

13.18 

1.15 

1.52 
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X  SQRT(X  +  5)  ♦  1 

LIMIT  . 

X-*+oo  3 

SQRT( 4  X  ♦  1)  ♦  X 


1 

$QRT(4) 


TAN(X) 

LIMIT  . - 

X--J"  LOG( COS(X ) ) 


Ml  NF 


(SQRTC 2  X)  -  2) 

LIMIT  . 

X— 2  LOG(X  -  1) 


1.93 


3.63 


3.83 
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TABLE  I  I 

PERFORMANCE  OF  WANDERER 


Integral  computed 


Resul t 

equivalent  to 


00  2 


+  A  X  ♦  B 


2 

+  10  X  +9 


+  A  X  ♦  B 


2 

♦  10  X  +9 


dX 


dX 


0  (principal) 


L0GC2) 


*PI  B  ♦  3  %PI 
12 


3  A  L0GC3)  +  (B  +  . 3 ) % P I 
24 


2 

(X)  -  SIN(X)  dX 


%P  1 
%PI 

20  SJNUPI/20) 


Time  In 
Sec. 

0.55 

2.02 

7.42 

11.32 

1.67 

1.08 
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f: 

L 


3  D  -5  SQRT(X) 

X  %E  dX 


oo  2 

COS  (X) 


dX 


U  2 

X  +  49 


I 


dX 


SQRTCX)  (2  X  ♦  1) 


l 


o»  1/3 

X  dX 


(4  SQRT(X)  +  3) 


3/2 


D 

(2016  %E  )  /  78125 


SPI/28 


*PI/SQRT(2) 


DIVERGENT 


0.95 


1.39 


0.81 


6.62 
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APPENDIX  D 
NOTATIONS 

INF  plus  Infinity,  MINF  minus  Infinity, 

%\  V-l  *PI  X 

$E  base  of  natural  logarithm 

*  Multiplication  operator 

**  exponentiation  operator 

SQRT  square  root  operator 

PLOG(X)  %?\  S  lm(PL0G(X) )>  -%PI,  the  principal  branch 

LOG(X)  when  X  >  0  for  PLOG(X) 

GLOG(X)  2  %?\  >  lm(GL0G(X) )  0 

ZP  the  zero  parameter  (defined  In  chapter  4) 

BETA  the  beta  function 


GAMMA 


the  gamma  function 
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APPENDIX  E 

AN  UNDECIDABILITY  RESULT 

Let  SI  =  {P(X1,X2, . . .  ,Xn)}  be  a  set  of  polynomials  with 
Integral  coefficients  In  XI,  ...  ,  Xn  and  S2  a  set  of 
functions  F  of  the  form 

F(X1,  ...  ,  Xn) 

2  2  tv  2  2 

=  (n  +  1)  (P  (XI,  ...  ,  Xn)  +  X  (SIN  *Xj)KJ  (XI,  ...  ,  Xn)>-1, 

where  Kj  Is  the  dominating  function  |6|  for 

d  2 
—  P. 

dXj 

Richardson  has  shown  |22,6|  that 

Lemma  1.  For  F  In  S2  the  predicate  "there  exists  an  n-tuple 
B  of  real  numbers  such  that  F(B)  <  0"  Is  recursively 
undec I dable. 

Lemma  2.  If  F(B)  <  0  for  some  n-tuple  B  real  numbers  then 
there  exists  an  n-tuple  A  of  nonnegative  Integers  such  that 
P( A)  =0  and  therefore  F(A)  =  -1. 
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Lemma  3.  Let 

3 

h(X )  =  X  SIN(X)  and  g(X)  =  X  SIN(X  ) 

Then  for  any  real  numbers  Al,  ...  ,  An  and  any  0  <  6  <  1 
there  exists  b  >  0  such  that 

|h(b)-Al|<6/  |h(g(b))-A2|<6,  ,..  ,  |h(g(...g(b)...)) -An  1  <€, 

By  use  of  these  Lemmas  we  can  show 

Corollary.  For  any  G  In  a  set  of  functions  of  the  form 
G(X)  *  F(h(X)/h(g(X))#  ...  ,  h(g(  • . .  (g(X) ) . . . ) > )  +  1/2, 
the  predicate  "there  exists  a  real  number  t  such  that  G(t)  a 
0"  Is  recursively  undecldable. 

Proof;  Suppose  the  predicate  Is  recursively  decidable/  then 
we  have 

1)  If  there  exists  a  t  such  that  G(t)°0  then  there  exist 
real  numbers  Al/  ...  ,  An  such  that  F(A1/  ...  ,  An)  <  0; 

2)  If  there  exist  real  numbers  Al/  ...  /An  such  that 
F(A1/  ...  /  An)  <  0  then  (by  lemma  2)  there  exist 
nonnegative  Integers  Bl/  ...  /  Bn  such  that 

F(B1/  ...  /  Bn)=-1.  Thus,  there  ex/sts  a  real  number  c  such 
that  G(c)  <  0.  Since  G(X)  can  be  large  and  positive  and 
GCX)  Is  continuous/  this  Implies  that  there  exist  a  t  such 
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that  G(t)»0. 

This  contradicts  lemma  1. 

From  the  definition  of  the  Function  F,  one  can  see  G  is 
always  -1/2  and  G  Is  large  and  positive  except  at  the 
vicinity  of  a  finite  number  of  points  where  G  Is  negative. 

Theorem;  The  convergence  of  a  set  of  Integrals  of  the  form 

dX 

2  2 
+1)G  (X) 

Is  recursively  unde'; Idable. 

Proof:  This  Integral  Is  convergent  If  and  only  If  G(X)  has 


no  real  zero 
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